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CHAPTER 1 


1.1 


(a) P 


P 

RT 


19 x 10 4 
(287)(203) 


0.326 kg/m ¿ 


(b) T 


P 

pR 


1058 

(123 x 10 -3 )(1716) 


501 °R 


1.2 N' = - | (p u eos 9 + t u sin 0) ds u 
r w 

+ J (p f eos# — T f sm#) ds f 

ds eos 9 = dx 
ds sin 0 = -dy 

Henee, 

N'= - (Pu-P{)dx + (x u + x f ) dy 

N' = - [(Pu-P«)-(P<-p»)]dx+ (r u + x ( ) dy 


Divide by q„ S = q„ c(l) 


N’ 1 rTE 

q*c c '“ 



v q„ J 


„ y 
Pf-P. 


V q m 7 



TE 

LE 




This is Eq. (1.15). 


(1.7) 
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r TE 

A' = I (-p u sin© + t u cosG) ds, 

J LE 


+ I LE ^ ? + T C C0S ®) ( (1 - 8 ) 

A ' = (Pu-p e )dy+ (t„ + z e ) dx 

A ' = J LE [(Pu-p=o)-(p f -p«)]dy+ |J (Tu + T f )dx 


A’ _ 1 p TE 

q„c c-’le 


Pu ~P. 

v. q* 


r 


p t -p. 

v q« 


dy + i f‘íiL + A 

c. ° V n n 


dx 


= “ íZ (c p,- c p ; ) d y + 7 r K-c«)dx 


Thisis Eq. (1.16). 

r te 

M'le ~ I [(p u COS0 + x u sin9)x - (p u sin© - x u cos0)y] ds u 

-í l£ 

+ f [-p, eos© +1 1 sin0)x + (p É sin© + x f cos0)y] ds ¡ 


M' le = [p„-p ( ]xdx- (xu + x^xdy 

+ tPu-p f ]ydy+ (x« + x f ) y dx 

r te r TE , . 

M'le= I [CPu * P») - (p í -Poo)]xdx- I (x u + x f )xdy 

J LE ^ LE 

+ 17! EP»-P~)-(P ( -P»)] y# + Í. P (x u + x f )ydx 

■/ LE w Le- 


Divide by q^c 2 : 
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q.c- 



X 

Pe -P» 

V q„ ) 



x dy 
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M’ 


^Cp ' 


LE 


-(Pí“pJ- 


N' 


(Pí-Pu) C 


( Xcp = c/ 2| 


1 For a flat píate, 0 = 0 in Eqs. (í .7) - (1.11). Henee, 

N' = P (p f -pu) dx= P [-2 x 10 4 (x-1) 2 + 1.19 x 10 5 ] dx 

J O J Q 


N r = -2xl0 4 [y-x 2 + x]i + [1.19 x 10 5 x] ’ = 1 i.12x10^N 


A' = f C (x r x u )dx= í' (731 x "~°' 2 + 288 x -0 ' 2 ) dx 

Jo Jo 


A' = [1274 x° £ } ’ = 11274 N 


L/ =N f cosa - A' sina = 1.12 X 10 5 eos 10° - 1274 sin 10° 


1.105 x ÍO-'N 


D' =N' sina + A r cosa= 1.12 x 10^ sin 10°+ 1274 cosa 


= £.07x10 N 


M' LE = Í C [p u -p f ]xdx = P [2 x 10 4 (x-1) 2 —1.19 x 10 5 J x dx 

Jo J O 


+ 2xl0 4 [—- — + [0.595x10^1 == 1-5.78 xl0 4 Ñm 

4 3 2 


M'c /4 = M' le + L' (c/4) = -5.78 x 10 4 + 1,105 x 10 5 (0.25) 
= -3.02 x 10* N/m 


M’je 

x cp =-— = - 


N’ 


(-5.78 x 10 4 ) 
1.12 x 10 5 


0.516 rd 
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1.5 


c = Cn cosa - c a sma 


: (1.2) eos 12° - (0.3) sina = |L18j 


C(j = c n sina + c a cosa 


= (1.2) sin 12° + (0.3) cosa = ¡0.279 


1.6 c n = c ( cosa + Cd sina 


Also, using the more accurale N' rather than L' in Eq. (1.22), we have 


x 


cp 


c M' c/4 
4 N’ 


c 

4 


r 

c 

V 





J 


Henee; 


a(°) 

c n 

Xcp/c 

-2.0 

0.0498 

1.09 

0 

0.25 

0.41 

2.0 

0.44 

0.336 

4.0 

0.639 

0.306 

6.0 

0.846 

0.293 

8.0 

1.07 

0.284 

10.0 

1.243 

0.277 

12.0 

1.402 

0.271 

14.0 

1.52 

0.266 
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Note that Xcp moves forward as a is increased, and that it closely approaches the quarter- 
chord point in the range of a of 10° to 14°. At higher angles-of-attack, beyond the stall (a > 
16°), x cp will reverse Its movement and move rearwaxd as a continúes to increase. Compare 
the above variation with the center-of-pressure measurements of the Wright Brothers on one 
of their airfoils, shown in Fig. 1.28. 


1.7 K = 3 (mass, length, and time) 

fi (D. pt,. V», c, g) = 0 Henee N = 5 

We can write this expression in terms of N — K = 5 — 3 = 2 dimensíonless Pi producís: 

Í2 (Hi, II 2 ) 


where 


fli = f 3 .(p^ c, D) 
n 2 ti (peo; Voo, c, g) 
Let n, = p OT a Vj c d D 
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l = (m r 3 ) a (/ f’) b £ c (m£ f 2 ) = O 
mass: a + 1 = O 
length: -3a + b + c + 1 = O 
time: -b - 2 = O 


Henee: 


D 


n x -orn T : 


D 


A.V.V 


1 


/>.Y.V 


D 


n,= 

q. ci 

Let n 2 = p ra 2 V K c b g d 

1 = (m r 3 ) a (£ f 1 ) £ b (¿ t' 2 ) d = O 

mass: a=0 

iength: -3a +l+b+d = 0 
time: -l-2d = 0 

Henee: 


n 2 = 


v, 

-fig 


Thus: 


f 2 (n Ií n 2 )=f 2 


f d v A 


U « c2 ’yfcg. 


= 0 


on 


|Cp = f(F r )| 


a = -1 

b = -2 

c = -2 


a= O 


d = -1/2 
b = -1/2 


1-8 Dw fl CpcOa C, Cp, Cy) 

K = 4 (mass, length, tíme, degrees) 
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í¿ (Dw> p* 5 V», c, a», Cp, c v ) O 

Henee, N — 7. This can be written as a fimetion of N — K = 7 — 4 = 3 pi producís: 

Í3 = ( 111 , IT 2 ,II 3 ) = 0 

where: 


Oí £1 (peo, Vce*? e, Cp, D) 


n 2 = fs (p ro , Veo, C, Cp, a*) 
n 3 = U (peo. Veo, c, Cp, Cv) 


The dímensions of Cp and Cv are 


[Cp] = 


energy 

mass(°) 


(force)(distance) (mft~ 2 )(■£) 
mass(°) m(°) 


[c p ] = i 2 f 2 O* 1 where (°) degrees. 


For np 

P JvJc k Cp n D = n, 

(m r 3 y (£ f’y (tf (£ 2 t' 2 T (V ( m £ t_2 ) = 1 


mass: i + 1 ~ 0 i -1 

length: -3i + j +k + 2n + 1 = 0 n=0 

time: -j — 2n —.2 ”0 j 2 

degrees: -n = 0 k~~2 

Henee: 


n,= 


D 


c- 


orlli 


P 

q«,c 


For n 2 : 

nh^pJ V w (7 c p k a* n 
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i = (m r 3 y {i f 1 ) nin 2 t 2 f (°y k (j r 1 ) n 

i = O 
k = O 
n = -l 
j = 0 


i = (m r 0 )' (i t - l y e (£ 2 t- 2 y (V a r 2 ) (°y’ 

mass: i = O 

lengtk -3i+j+k + 2n + 2 = 0 
time: -j - 2 n - 2 = O 

degrees: -n - 1 = O 

Henee: 

c 

II 3 = ~. We can take the reciprocal, and stíll have a dimensionless product. 

c p 

Henee, 

n 3 =^=y 

C P 

Thus, 

í D V C n ) 

£ ^ co p 

r 3 o 3 » 

vq«,e a „ c v 2 



mass: i = 0 

length: -3i + 1 +j + 2k + n = 0 
time: -1 - 21 c - n = 0 

degrees: -k = 0 -2 

Henee: 



For n 3 : 

n 3 — ppo' Vj c k c p n c v 
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lC D = f (M^, y H 


19 E = M PRj 

M 2 V 2 a, V^T, 200 V 200 
Henee, the Mach numbers of the two flows are the same. 

Re, p,Y,c, f a a p.v.c, ItT rmViooyn |800 = , 0 3 54 

R e¡ p V,c 2 lp,j p V.c. V T, O.739A200A2Jll200 

The Reynold’s numbers are difieren! . Henee, the two flows are not dynamícally similar. 


L10 Denote free flight by subscript 1 ? and the wiud tunnel by subscript 2. For the lift and 
drag coefficients to be the same in bolh cases, the flóws must be dynamically similar . Henee 

Mj — Mo 


and 


Re] — Re2 


For Mach number: 



a 


1 


a 


2 


Since a a a/t , we have 


V 2 _ V, 25 Q 

vt'Vt "V223 


-16.7 


( 1 ) 


For Reynolds number: 


. P. v i c , 


Mi 


¿b V 2 C 

Mi 


Assume, as before, that fr a Vt . Henee 


or 5 

P 2 V 2 = £&_(O = (0.414X250) ^5^1 
^ W 223 U 

' or. 


P2Y2 


= 34.65 


Finally, from the equatíon of State: 


__ p 2 1.01 x 10 
p? 1 2 --- 

R 287 


: 351.9 


( 2 ) 


(3) 


Eqs. (1) - (3) represen! three equations for the three uhkno-wns. p 2 , V 2 , and X 2 - They are 
summarized below: 


V, 


= 1.67 


0 ) 


P2X1 


= 34.65 


( 2 ) 


P 2 T 2 = 351.9 (3) 

FroraEq. (3): 

p2 = 351.9/T 2 (4) 

Subst. (4) into (2): 


351.9 í V 2 

t 2 jt) 


= 34.65 (5) 


Subst. (1) into (5); 


3519 

mí- (16.7) = 34.65 

F. 


Henee, 


n 


t _ (351.9)(16-7) _ 
T2 “ (34.65) “ 


169.6°K 


From Eq. (1): V 2 = 16.7 ^ = 16.7 Vi 69.6 = 


217.5 


m 


sec 


From Eq. (3): p 2 : 


351.9 351.9 


169.6 


2.07 


kg 


1.11 Pb = Pa-pgAh 


1.01 x 10 5 ~ (1.36 x 10 4 )(9.8)(0.2) 


p b — 17.43 x 10 4 N/m 21 


1.12 Weight = Buoyancy forcé + lift 


w « 

B + L 

B = (15,000) 

(1.1117) 

(9.8) = 1.634 x 10^ 

V 

acceleratíon 

volume 

air density 

(m 3 ) 

at 1 OOOin 

of gravjty 


(kg/m 3 ) 

(m/sec) 

__ 1 2 
q» Peo Voo 

2 

= - (1.1117) 

2 

(3 O) 2 = 500 N/m 2 


S = rai 2 /4 = tc( 14) 2 /4 = 153.9 m 2 
L = q«, S C L = (500)(153.9)(0.05) = 3487 K 


Henee: 


W = 1.634 x 10 5 + 3847 = |l.67x 10V4 


1.13 Let us use the fonnalism surrounding Eq. (1-16) in the text. In this case, c<¡ c a , and 
firomEq. (1.16), neglecting skin friction 
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1 

Cd= - 
c 


( 1 ) 


J™ K-c,> 


From Eq. (L13) in the text Eq. (1) above can be written as 

Cd= 7 C ( C p, ~ C p,) (“ sin0ds ) O) 

Draw a picture: | 

-© i 

Following our sign 
convention, note 

that 0 is drawn V¿o 

counterclockwise - 

in this sketch, henee 
it is a negative 
angle, -0. 

From the geometry: 

- 0 = 7t - <f> 



Henee, sin (-0) = - sin 9 = sin (7i-0) = eos <j> 

Substitute this into Eq. (2), noting atso that ds = rd<¡> and the choxd c is twice the radius, c = 
2r. From Eq. (2), 

Cd= ¿ ^ (C p „ — C Pc ) eos <j> r d <|> 

C ( C p. _C p t ) cos, t>d<¡> 

Cd= Í C Po eos <(> d <|i ~ C p¡ eos <t> d <j> (3) 


Consider the limits of integration for the above integráis. The first integral is evaluated from 
the leading edge to the trailing edge along the npper suiface. Henee, <j> = 0 at LE and % at TE. 
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The second integral is evaluated from the leading edge to the trailing edge along the bottom 
surface. Henee, <j> = 2n at LE and k at the TE. Thus, Eq. (3) becomes 


Cd = 

1 [' C 

2 J o * 

1 r x 

eos 6 d <b — C 

> ? J 2* Pí 

In Eq. (4), 



C p„ 

= 2 eos 2 <j> 

for 0 < <j) < 7i/2 

C p, 

-0 

for y < 9 < ti 

C P t 

= 2 eos 2 <j> 

3 71 

for — <§<2tc 

C p, 

= 0 

_ , 3tT 

for 7 r < 9 < — 

Thus, Eq. (4) becomes 



f o 

r 3jz72 -5 

Cd = 

eos 

d 0 

<j> d <j> - eos (j) d <j> 


1 2 

Since eos <¡» d <¡> = (— sinocos <j> + 2), Eq. (5) becomes 
3 

c d = [(^- sin<j))(cos 2 (j) + 2)] f- - [(- sin<|))(cos 2 (j) + 2] 
J> i 

Cd=(^)(l)(2)-(|)(-1)(2) 

Cd - 4/3 


1.14 


(4) 
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Consider the arbitrary body sketched above. Consider also tbe vertical cylinder element 
inside the body which intercepts the surface area dAj near the top of the body, and dAz near 
the bottom of the body. The pressures on dAi and dA2 are p 3 and P2 respectively, and makes 
angles 0] and 02 respectively with respect to the vertical line through the middle of dAi and 
dA2. Tbe net pressüre forcé in the y-direction on this cylinder is: 

dF y = - pi eos ©i dAj + p2 eos 02 dA2 (1) 

Let dA y be the projection of dAt and dA2 on a plañe perpendicular to the y axis. 


dA y = eos ©i dAi = eos 02 dAi 


Thus, Eq. (1) becomes 


dFy=(p2-p))dA y 

From the hydiostatic equation 
P2 -Pi“ pgdy 

J “i 


( 2 ) 


( 3 ) 


Combining Eqs. (2) and (3), 

dFy^ í* 12 pgdydAy (4) 

J h l 

However, dy dA y = dV = element of volume of the body. Thus, the total forcé in the y 
direction. F y; is given by Eq. (4) integrated over the volume of the body 
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f, = 

'^r' —v- 

Forcé on body Weight of fluid displaced by body. 


1.15 From Eq. (3.45) 

_ L _ 2W 2(2950) 

CL ~ 9^S _ /^Vjs ' (0.002377)V o 2 (174) 


C L = 


14265 



Also, 


C D = 0.025 +0.054 C L 2 (2) 
Tabulate Eqs. (1) and (2) versus velocity. 





L C L 

(fl/sec) 

C L 

Cu 

d'c d 


70 

2.911 

0.483 

6.03 

90 

1.761 

0.192 

9.17 

110 

1.179 

0.100 

11.79 

130 

0.844 

0.063 

13.40 

150 

0.634 

0.047 

13.49 

170 

0.494 

0.038 

13.0 

190 

0.395 

0.033 

11.97 

210 

0.323 

0.031 

10.42 

230 

0.270 

0.029 

9.31 

250 

0.228 

0.028 

8.14 


These results are plotted on the next page. 
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5 


¿o- 

L 


o 

¡A- 

)4- 

/,?■ 

ti- 

/.O 

A5 - 

0,6- 

8" 

0,0 

6- 

0.4- 

f- 

0,2- 

z- 

O 

o - 


L ^V 



^~o 

Ó ./? 

a/2 


o. to 


0,0 s 

O, D O 

0.04 

0 . 0 % 


OÓ 


—I-:-f-1-1- 1 -1-•— 

90 /JO /SO /SO /7D /?0 

i/£í-oc/t y, l¿: 

J 


- i — 

¿?/o 


--,-L. o 

^3a 


C f^r/s^rc ) 


Examining tliís graph, we note, for steady, level flight: 


1. The lift coefficient decreases as V m íncreases. 

2. At lower velocity range, the drag coefficient decreases even faster fhan the lift 
coefficient with velocity. (Note that on the graph the scale for Cd is one-tenth 
that for Cl-) 

3. As a result, the lift-to-drag ratio first increases, goes throngh a máximum, and 
then gradually decreases as velocity increases. 


It can be shown that the máximum velocity for this airplane is ábout 265 ñ/sec at sea level. 
As seen in the graph, the máximum valué of L ÍD occurs around V<*> — 140 ft/sec, which is 
much lower than the máximum velocity. However, at higher velocity the valué of L/D 
decreases only gradually as V„ increases. This has the practícal implication that at higher 
speeds, even though the valué of L/D is less than its máximum, it is stíll a reasonably high 
valué. The range of the aircraft is proportional to L/D (see for example, Anderson, Aircraft 
Performance and Pesien. McGraw-Hill, 1999, or Anderson, Introduction to Flight 4* ed.. 
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McGraw-HilL, 2000), To obtain máximum range, the airplane should fly at the velocity for 
máximum L/D ? which for this case is 140 fl/sec. However, one reason to fly in an airplane is 
to get from one place to another in a reasonably short time, By flying at the lov¿ velocity of 
V* - 140 ft/sec s the flíght time may be unacceptably long, By cruising at a higher speed, say 
200 ft/secx tb.e flight time will be cut by 30% ? wíth only an 18% decrease in L/IX 
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CHAPTER 2 



F = -^pdS 


If p = constant = p» 
F'-P.^páS (1) 


However, the integral of the surface vector over a closed surface is zero, i.e., 

^di = 0 

Henee, combmmg Eqs. (1) and (2), we have 


F - 0 


2.2 



/~^7 7 7 7 7 / / 7 7 

¿-owei- Wa/¡ fa. 


*L- 
7 ~y 


19 





Denote the pressure distnbutions on the upper and lower walls by p u (x) and p t (x) respectively. 
The walls are cióse enough to the model such that p u and p ( are not necessarily equal to p^ 
Assume that faces ai and bh are far enough upstream and downstreara of the model such that 

p — p« and v = 0 and ai and bh. 

Take the y-component of Eq. (2.66) 

L = -§ (/?V ‘ dS) v - J{ (pdS)y 

S abhi 

The first integral = 0 over all surfaces, either because Y' ds = 0 or because v = 0. Henee 

b h 

[ J Pu dx - j p f dx] 

Minus sign because y-component is in downward 
Directíon. 

Note: In the above, the integráis over ia and bh cancel because p = p TO on both faces. Henee 


h 

b 

L' = f p £ dx- 

í Pudx 

i a 


L’ = - JJ (p¿S)y =- 

ábht 


dy _ v __ cy/(x 2 +y 2 ) y 
dx u cx/(x 2 +y 2 ) x 

dy _ dx 
y x 

£n y = f n x + Ci = f n (C 2 x) 


ly^czxl 

The streamlines are straight lines emanaíing from the origin. (This is the velocity field and 
streamline pattem for a source, to be discussed in Chapter 3.) 


94 dy _ v _ x 
dx u y 


y dy = - x dx 
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y 2 = 'X 2 + const 
x 2 + y 2 = const. 

The streamlines are concentric with their centers at the origin. (Thís is the velocity fíeld and 
streamlíne pattem for a vortex. to be discnssed in Chapter 3.) 


2.5 From inspectioiij since there is no radial component of velocity, the streamlines must be 
circular, with centers at the origin. To show this more precisely, 

u = - Ye sin = - cr — = - cy 
r 

v = V e eos 9 - cr — = ex 
r 

dy _ v _ x 
dx u y 


y r + x 2 = const 

This is the equation of a circle with the center at tire origin. (This velocity field corresponds to 
solid body rotation.) 


2.6 *L = 1=-1 

dx n x 

dy dx 
y x 

(ny = x (nx + Cj 
y = c 2 /x 

The streamlines are hipérbolas. 
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2.7 


(a) 


<5v Dt 


In polar coordinates: 
Transformation: 


V ‘^ = rl (rVr)+ r^ 


x = r eos 0 


y = r sin 0 

V r = u eos 0 + v sin 9 
Ve = - u sin 9 + v eos 0 
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ex 


u = 


cr eos 9 c eos 0 


(x 2 +y 2 ) 


v = 


cy 


cr sin# c sin# 


(x 2 + y 2 ) 


V r = — cos 2 0+ - sin 2 0 = — 

r r r 


c c 

Ye = — cos0 sin9 + — cos9 $in9 — 0 
r r 


y 24 «0 + 132»-o 

r de r 09 


(b) From Eq. (2.23) 


VxV=e í 


& r r dQ 


VxV = e z [0 + 0-0] = 
The flowíield ís irrotational. 


2.8 u= Cy 


cr sin# c sin# 


(x 2 +y 2 ) 


v = 


-ex cr eos# c eos# 


(x 2 +y 2 ) 


c c 

Y r — co.sQ sín0 — cos9 sin@ = 0 
r r 


Ve = — sin 0-cos^0 =- 

r r r 


(a) v-V = Í 4 ( 0 ) + I ^£)-0 + 0-§ 

r ¿k r ov 



(b) V x V = e 2 


¿?(-c/r) c 1 ¿7(0) 

<3: r 2 r 


c c 

= e z —-0 


Vx V=Ó¡ except at the origín, where r = 0. The flowfield is singular at the origin. 




Mass flow between streamlines ~ Ay/ 


Ay/ =p V An 


Ay/ - (-p V 0 ) Ar + p V r (r0) 
Let cd approach ab 

dy = - p Ve dr + p r V r d0 

Also, since yz = y/ (r,9), from calculus 


d v ,AL <k+ íy/ 

g¡ se 


Comparing Eqs. (1) and (2) 


de 


-pV e = 


dy/ 


and 


prY r = 


dy/ 

He 


or: 


pV r = 


1 dy/ 
r dO 


P V e = - 


dy/ 

IT 


2.11 u = ex = = cxy + f(x) 

ds 


v = - cy = - : y = cxy + ffy) 

dn 


( 1 ) 


( 2 ) 


( 1 ) 

( 2 ) 
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Comparing Eqs. (1) and (2), f(x) and f(y) = constant 


|v|/ = c x y + const. | 

n = ex = : <j> = ex 2 + f(y) 

dx 

v = - cy — : <j> = - cy 2 + f(x) 

Comparing Eqs. (4) and (5), f(y) = - cy 2 and f(x) = ex 2 
4> = c (x 2 - y 2 ) 

Drfferéntiatmg Eq. (3) with respect to x, holding \\i — const. 

n dy 
0 = ex — + cy 
dx 

or, 

f) 

y =const 

Differentiating Eq. (6) with respect to x, holding <j> = const. 

dv 

0”2cx-2cy— 
dx 


Oly 



-x/y 

0=con sí 


Comparing Eqs. (7) and (8), we see that 

M =_ 1 _ 

(¿i) 


Henee, lines of constant are perpendicular to lines of constant <j>. 




2.12. The geometry of the pipe is shown below. 



* 


As the flow goes through the U-shape bend and is tumed, it exerts a net forcé R on the infernal 
surface of the pipe. From the symmetric geometry, R is in the horizontal direction, as shown, 
acting to the right. The equal and opposite forcé, -R, exerted by the pipe on the flow is the 
mechanism that reverses the flow velocity. The cross-sectional area of the pipe inlet is rrd 2 /4 
where d is the inside pipe diameter. Henee, A = 7td 2 /4 = -7t(0.5) 2 /4 = 0.196m 2 . The mass flow 
entering the pipe is 


m =pi A Vi = (1.23)(0.196)(100) = 24.11 kg/sec. 

Applying the momentum equation, Eq. (2.64) to this geometry, we obtain a result similar to Eq. 
(2.75), namely 


R = -§(pVdS)V (1) 

Where the pressure term in Eq. (2.75) is zero because the pressure at the inlet and exit are the 
same valúes. In Eq. (1), the product (p Y ' dS) is negatíve at the inlet (V and dS are in opposite 

directions), and is positive at the exit (V and dS) are in the same direction). The magnitude of p 

♦ 

Y ' dS is simply the mass flow, m. Finally, at the inlet Vj is to the right, henee it is in the 
positive x-direction. At the exit, V 2 is to the left, henee it is in the negative x-direction. Thus, 
V 2 = _ V[ With this, Eq. (1) is written as 


R = - [-m V! + m V 2 ] = m (Vi — V 2 ) 


= m[V 1 -(-V,)] = m(2V 1 ) 
R= (24.11)(2)(100) = [4822 N 
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CHAPTER 3 


3.1 Consider steady, inviscid flow. 


x-moxnentum: 


y-momentum: 


di di di da 

pu — ■+■ pv — + pw - 

P da dy dz da 


di di di dp 

pu — + P v — + pw — - - —- 

da dy dz dy 


( 1 ) 

( 2 ) 


dx dx dx _ dp 

z-momentum: pu + pv + pw —-— 

d& dy dz dz 


(3) 


Multiply (1), (2), and (3) by dx, dy, and dz respectively: 


¿b di , di , 1 dp , 

u—dx + v — dx + w — dx-dx 

da dy dz p da 


(4) 


di , dv di _ 1 dp 

u — dy + v — dy + w — dy-—- dy 

da dy dz p dy 


dx , dx dx __ I dp , 

u — dz + v — dz + w dz - - ■-— dz 

da dy dz p dy 


(5) 

( 6 ) 


Add (4) + (5) + (6): 


(di dv dx A (di, di, dx 
u —dx+ — dy +—dz +v —-dx + ~dy + —-dz 
da da ) \dy dy dy J 


( di , di , <5W , ^| 1 

+ w —dx +—dy+ —dz =- 

\dz dz dz ) p 


— dx + — dy + -^dz 
(da dy dz 


(J) 


For irrotational flow (see Eq. (2.119)): V x V — 0 
Henee: 


dx di di _dx di _ di 
dy dz’ dz da’ di dy 


( 8 ) 


Subt Eqs. (8) into (7): 


2S 


u 


3a _ ¿Si , ¿Si , 

-— dx -i-dy + — dz 

K3í 3$ dz 


^ /* 
+ v 


j 


3v . ¿ 7 / . ¿ 7 / , 

— dx -h —dv + — dz 
^¿3 í dy dz ) 


, ¿Sv , ¿5W , ¿Sv _ ^ 1 

+ w I-dx -i-dv h -dz ( = 

dz dz ^ dz 


P 


^dX^dy + ^dz] 
¿fe dy dz. ) 


u du + v dv + w dw = - — dp 

P 


— d (u 2 + v 2 + w 2 ) = — d (V 2 ) = VdV = - — dp 

2 2 p 


dp = - pV dV which mtegrates to 


p + ■— p V 2 = const. 


for incompressible flów. 


3.2 V, = 


(2)(p, -p 2 ) 


P\ 


ívV 

\A,, 


-1 


pj = 2116 lb/ft 2 ; p2 = 2100 lb/ft 2 , A 2 /Ai = 0.8 


V,= 


2(2116 

- 2100) 

(0.002377) 

[m'-il 

vO.S/ 


154.7 ft/seq 


3.3 pi~P2= Vfe p V, 2 


A, 
A 2 J 


-1 


= Vz (1.23)(90) 2 [(1/0.85) 2 - 11 =|1913 N/m" 


3.4 


V - 2 w A h 

y f\} — (A 2 / a| ) 2 ] 
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w= p m g = (1.36x 10 4 ) (9.8 —^-) = 1.33 x 10 5 N/m 2 

sec 


Ah = 10 crn = O.lm; p = 1.23 kg/m 3 , 


Ai 

A, 


12 


2(1.33 

x 10 f )(0.1) 

(1.23) 




1,147 m/sed 


3.5 Pi ~ Pi = w Ah = (1.33 x 10 5 )(0.1) = 1.33 x 10 4 N/m 2 

p 2 = p, - 1.33 x 10 4 = 1.01 x 10 5 - 1.33 x 10 4 - 8.77 x 10 4 N/m 2 

p 0 = P 2 + ~ P V 2 2 = 8.77 x 10 4 + ~ (1.23)(147) 2 = 

Note : It mákes sense that the total pressure in the test section would equal one atmosphere, 
because the flow in the tunnel is drawn directly from the open anibient surroirodings, and for 
an inviscíd flow, we have no losses between the inlet and the test section. 


1.01 x 10" N/m 2 


3.6 


1 

Po = p® + — 

2 


pV ro 2 




|2(p„-p ffl ) _ ] 

Í2(L07 —1.01) x 10 5 

98 8 _ m 

P v 

123 

sec 


3.7 


C p = l- 


JV 


= 1 - 


130 

98.8. 


= 1-'0.731 


3.8 


< 

s 

-^ 
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It is a physícally possible incompressible flow. 


V x V 


= i (0-0)- j (0-g) + k (0-^) 


Vx V = 0 


i j k 

<3c djf dz. 

0 0 
u ^ y( 

The flow is irrotational. 


3.9 For a source flow. 


A 

y ^ y r e r = — e. 


2;zr 


In polar coordinates: 


r¿h r <5P 




r r dQ 



Henee, the flow is a physical possible incompressible flow, except at the origín where r = 0. 

What happens at the origin? Visualize a cylinder of 
radius r wrapped axonnd the líne source per unit depth 
perpendicular to the page. The volume flow across this 
cylindrical surface is 
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(I) 


§ V - as 

s 

Since we are considering a unit depth. then we have the volume flow per unit depth. This is 
precisely the definition of source strength, A. Henee, from (1), 

A = constató = ^ V ' dS (2) 

s 

From the divergence theorem: 

§ V - dS = (V-V) dV (3) 

S v 

Combining Eqs. (2) and (3) 

<£ji[ (V-V) dF= A = constató (4) 


Shritót the volume to an infinitesimal valué, AV, around the origin. Eq. (4) becomes 

(V-V) AV= A 

Taking the limit as AH-» 0 

(V - V) = lim —= oo. 

V ' AV 


Henee V-V = oo at origin 


To show that the flow is irrotational, calcúlate V *V. 



32 



Henee, 


V xV = O eveiywhere. 


.10 


<(> = Vocx; 


^=V„- 

o 1 <j> 

di 

dy 1 

O 

fl 

*1 

d 2 <j> 

dy 

dy 2 


= 0 


Henee, Laplaces equation: 


¿V , <? 2 


2 + ~ 0 + 0 = 0 is identically satisfied. 


ac ay 1 


Similarly, for y = V y: = 0, = 0 

di di 


&f ' ¿\' 2 


Henee, Laplaces eqnation: 


g! ^2 

——- + —= 0 + 0 — 0 is identically satisfied. 
d¿ dy 1 


3.11 <)> 

A 

Ai 

d<f> _ 

A 

1 d 


2 71 


’ a 

2 n 

r 




d(f> 

= 0, 

d 2 <j> 

= 0 




d0 


de 1 


Henee, Laplace 

’s equation 




i 

<?í r 

d$\ 

1 d 1 

K 

1 d 

[Xl 

r 

<3: l 

J 

r 2 dS 


2^1 

_OK _ 


+ 0 = 0 


is identically satisfied. 




<3- ¿i: 2 

de 2rv ’ d<l> 2 


Henee, Laplaces equation 


-l£í r A.') + -L£Í 

rál r 2 ¿39 


-^( 0 ) + ~( 0)=0 

r ¿i r 


is identically satisfied. 


3.12 The stagnation point is a distance A/27iV w upstream of the source. Henee, 
A 


2jrV 


= 1 , or A = 2n:V« 


The shape of the body is given by 


,, - A . A 

w - V^r sin 0 + — 0- 

2 % 2 


or. 


, . A A 
r sin 0 + -—— 6 = - 


or. 


2 tt\ 2V 


. „ 2ttV„ 2ttV„ 
r sin 0 +- 0 = - 


2r¿V 2V, 


or. 


Ir sm 0 T 9 = 7 i ¡ Equation of the semi-infinite body. 


r = 


k-B 
sin 9 
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9(rad) 


r 


x - r cos8 


71 

1 

3 

1.0033 

2.8 

1.02 

2.5 

1.072 

2.9 

1.255 

tt/2 

1.57 

1.3 

1.91 

1.0 

2.54 

0.75 

3.509 

0.5 

5.51 


To ülot the pressure coefficient ; 


y — r sinO 


-1 

0 

-0.990 

0.1416 

-0.961 

0.3416 

-0.859 

0.6416 

-0.522 

1.142 

0 

1.57 

0.511 

1.84 

1.372 

2.14 

2.57 

2.39 

4.84 

2.64 


Cartesian Coordinates of Body 


A 2 tN 

V r = Va,eos 9 + —— =V«cos0+—-= 
2ot 2^r 


= V„ eos 9 + 


V 


V e = -V w sin 0 


or. 


V 1 

—- = eos 9 + - 
V. r 


V, 


B - 


v_ 


- sin 9 


Í V ] 

'-(i'l 

2 

1 + 

ÍV ) 
v e 

IvJ 

IvJ 


{vJ 


2 i 2 

= eos (9 +— eos 6’ + — + sin 2 (9= 1 +—eos (9 + - 
r r r 











At point A: Velocity due to freestream = V» 


Velocity due to source = - 

2ír(r + b) 

(note that it is in the negative x-direction) 


, , , (+A) 

Velocity due to smk = —— - 

2tt(i + b) 

(Note that it is in the positive x-dírectíon) 


Total velocity at Point A: 


, r A 1 Al 

Ya = Ve__ -1___ 

2k (r - b) 2k (r + b) 


From point A to be a stagnation point. Va = 0. 


„ A r 1 1 , 

0 = V» + — [-+-]. 

2 % (r + b) (r - b) 


0 = V« 


A j ~ r — b - (r + b) = v A (~ 2b ) 


L (r + b)(r - b) 


2 tc r 2 - b 2 



3.14 


y = ^L 
di r 99 


_ , ,, k sin# 

For a doubletMj/ =- 

2 k r 



Substitute (2) into (1) 

d<¡) \f k coséR k eos 0 

di rl 2 n r / 2n r 2 

Integrating with respect to r 



or. 


, K C0S<9 

- - 

2 ti r 


2 \ 


3.15 \|/ = (V«rsmG)fl— 

V r ~ J 


V,= -^-=(V„c os0) 
r 39 


R 2>v 
V r 2 J 


"-f- 


1+4- 

V r 2 y 


V» sín9 


r r 2 V 

V 2 = V 2 + V e 2 = 1-Ar 

V r ' J 


V ro 2 cos 2 0 + 


r ! ' ! 

V l J 


V« sin 2 0 


V 2 ( R 2 V 

C p =l- —=1- 1-^- 

V.- \ T ) 

At the snrface, r = R 


cos 2 0 - í] _ R_^| 2 sin 2 0 


C p = 1 - 4 síxF0 


3.16 From Eq. (3.93): ^-= ^1 - ^ ^ 


l\ 


COS0 
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From Eq. (3.94): 


44 

V rV 


sin© 


At any given point (r,0), V r and Vg are both directly proportional to V» Henee, the dixection 

—^ 

of the resultante V 5 is the same, no matter what the valué of V„ may be. Thus, the s hape of 
the streamiines reniains the same. 





3.17 From Eq. (3.119): 



f T> - 'N 

1-4, 

v r; 


COS0 


From Eq. (3.94): 





i 

1 + — 
r j 


sin0 - 


r 

2*rV„ 


Note that VqÍVo, is itself a fimetíon of V«> via the second term. Henee, as V„ changes, the 
direction of the resultan! velocity at a given point will also change. The shape of the 
streamiines changes when V» changes. 


3.18 L' = p»V„r 

L' 6 


T- 


(1.23)(30) 


0.163 nT/sec 
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3.19 At standard sea level conditions. 


p M = 0.002377 

ft* 


M<* = 3.737x10' 7 


slug 

(ft)(sec) 


Also: 


V = 120 mph = 120 



ft/sec= 176 — 
sec 


q» =|p w V„ z =^ (0.002377) (176) 2 = 36.8 lb/ft 2 
For the struts: D = 2 in = 0. 167 ft. 


p _ pVT> _ (0.002377)(187.7)(0.167) 
p ~ 3.737 x 10 -7 


199,382 


From Fig. 3.39, C D = 1. The total frontal surface area of the struts is (25) (0.167) = 4.175 ft 2 . 
Henee, 


Drag due to struts: 

D s = q«S C 0 “ (36.8X4.175X1) = 153 Ib 

3 

For the bracing wires: D = — in = 0.0078 ft 

Re = 199382 í- Q ' Q - ° - 78 ] = 9312 
^ 0.167 ) 


From Fig. 3.39, Cd = 1- The total frontal surface area of the wires is (80) (0.0078) = 0.624 
ft 2 . Henee, 

Drag due to wires: 

Dw = q* S C D = (36.8)(0.624)(1) = 23 Ib 
Total drag due to struts and wires = D $ + Dw ~ 


153+23 = 176 
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The total zero-liñ drag for the airplane is (including struts and wires) 


. C Do = q« S C Do = (36.8)(230)(0.036) = ¡304.8 


Note thal for this example, the drag due to the struts and wires is 


176 

304.8 


= 0.58 of the total 


drag — í.e., 58 percent of the total drag . This clearly points out the drag reduction that was 
achieved in the early 1930’s when airplane designers started using intemally braced wings 
with one or more central spars, thus eliminating struts and wires completely. 


3.20 The flow over the airfoil in Figure 3.37 can be syntheized by a proper distribution of 
singularities, i.e., point sources and vórtices. The strength of the vórtices, added together, 
gives the total circulation, T, around the airfoil. This valué of F is the same along aü closed 
curves around the airfoil, even if the closed curve is drawn a very large distance away ffom 
the airfoil. In this case, the airfoil becomes a speck on the page, and the distributed point 
vórtices appear as one stronger point vortex with strength T. This is exactly equivalent to the 
single point vortex in Figure 3.27 for the circular cylinder, and the lift on the airfoil where 
the circulation is taken as the total T is the same as for a circular cylinder, namely Eq. 

(3-140), 


L' = V«>r 



r 


/ 






CHAPTER 4 


2 _ 1 tr\ f\rv^->n / zr\'\~ . 


4.1 q w = ~ p* V.* = ~ (0.002377)(50)" = 2.97 Ib/fr 


c f =0.64 and c m =-0.036 

* "V-í 


L' = q TO s c f = (2.97)(2)(1)(0.64) = [3.80 Ib per umt span 


M c/4 Qoo S C C m = (2.97X2)0)(2)(-0-36) = |-0.428 ft/lb per umt span[ 


^ - • 2_ 1 /-i ^'>\/-rr í \2 — iíoo\t/~2 


4.2 q. = - p» V» = ^ (1.23)(50) =1538 N/m 


L’ 


1353 


q w s (1538X2) 


0.44 


From Fig. 4.5, 


be = 21 




DV 

Dt 


ds + 


i 


V • ds 



Dt 


Henee, tbe second temí in Eq. (1) becomes 


r -+ j 

r ÍV 2S I 

cj V-dV = c| 

1 4 T 


From the momentum eqnation. 
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DV 1 

—- =-Vp (neglecting body forces) 

Dt p 


Henee, the fírst temí in Eq. (1) becomes 


4 


DV 

Dt 


"* r 1 r 

ds = -<j — Vp'ds=-<j 


dp 

P 


When p = const, or p = p(p), then 



(4) 


Substituting Eqs. (2) and (4) into (1), we obtain 

Note: See Karamcheti. Ideal-Finid Aerodynamics. for moré 
details (pp. 239-242). 



4.4 M'le = -p.V. f 

J O 


Ptc v» f' ~ (i - cose) e (y)^ sine de 

J ° 2 2 


-p 0 V»-2aV 0= [' (1 -cos 2 0) 

A J o 


de 


p.v»—a 
2 


_ , 1 V 2vj*« 

( 2 P» ) c 2 


= - q w c 2 — This is Eq. (4.36). 
2 


4.5 c, = 2 7t a where a is in radiaras. Henee 
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C, = 271 


1.5 

573. 


0.164 


c m = " Qf/r = I-Q.041 


4.6 (a) 


For0 < — <0.4: \—) =0.2-0.5 
c Vdx/, 


For 0.4 < — < 1: í—) = 0.0888 - 0.2222 [ - 
c vdxy, 


Since x = — (1 ~ cosG), then 


= -0.05 + 0.25 cos0, for 0 < 9 < 1.3694 


-0.0223 + 0.1111 cos0, for 1.3694 <0 <7i 


a w =-- Í' T — (cos0 - 1) d9 
K J 0 dx 


1 r 13694 1 r 

= — f (-0.05 + 0.25 cos9)(cos0 - l)d0-- f 

K J ° JC J 


1.3694 


(-0.0223 + 0.1111 cos0)(cos0 -1) d0 
| (0.05 - 0.3 cos0 + 0.25 cos 2 0) d0 

J o 


1 p 1.3694 
7T 


71 


f * (0,0223 - 0.13334 cos0 + 0.1111 cos 2 0) d9 

J 1.3694 7 


4 [0.059 - 03 sin.9 + 0.25 (- + - sin29) j 1 /® 4 
ti 2 4 
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4.7 


1 9 1 

- — [0.0223 9 - 0.1334 sinQ + 0.111 (- + - sin20)]^ 694 
n 2 4 

= -— [0.06847-0.2939 + 0.1712 + 0.0245]-— [0.0701 + 0.1745] 


+ 1 [0.0305-0.1307 + 0.0761 + 0.0109] 

K 


-0.2281 


=-0.0726 rad =4 4.161 


n 


(b) 

c f = 2 7i (cc +- ai=o) where a is in radíaos 


27 t 

Ü73 


[3-(-4.16)] = 10.782 


Aj = — f —cosG d0 
J 0 dx 


2 f 1.3694 

= — í (0.05 + 0.25 cose) cose 

7t J 0 


de 


+ —[* (-0.0223 + 0.1111 cose) 

K .> 1.3694 


cose de 


= — í 1 (-0.05 cos0 + 0.25 cos 2 0) de + 

7t J o 


— r * (-0.0223 cose + 0.1111 cos 2 e) de 

% j 13694 


- [-0.05 sine + 0.25 (- + - sm20)]’- 3694 
7t 2 4 


? 9 1 

— [(-0.0233) sinQ + 0.1111 (—* sin20))^ T 3694 
n 2 4 
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7T 


[-0.04899 + 0.25 (0.6847 + 0.09800) + 0.1745 


+ 0.02185 - 0.1111 (0.6847 + 0.09800)] 


A, = (0.2561) - = 0.1630 
% 


A 2 = - í" — cos20 d9 

TT ^ í* 


2 r 1.3694 2 r w 

F (-0.05 + 0.25 cosG) eos 20d0 + — I (-0.0223 

Jo 77 - J 1.3694 


71 


+ 0.1111 cosG) cosG d9 


— [- (-0.05) sin 20 + 0.25 (-55£ + 
n 2 2 6 


+ - (-0.0223) sin 29 + 0.1111 ( ^ + ^£)]* 

TT 2 2 6 


- ['0.009800 + 0.25 (0.4899 - 0.1372) + 0.004371 

71 


0.1111 (0.4899 - 0.1372)] 


(0.0436) - - 0.0277 
n 


71 


71 


c mc/4 = — (A 2 -Aj)= - (0.0277 - 0.1630) = [-0.1063 


| [1 + — (A, -A 2 )] = \ [1 + ^ (0.1630-0.0277)3 =g3g 


c 4 


4" 0.782 
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Experiment (Ref. 11) 

Theorv 

% Difference 

OCl=0 

-3.9° 

-4.16° 

6.25% 

C f 

0.76 

0.782 

2.8% 

C„ 

-0.095 

-0.1063 

10.6% 


4.9 M'le = - P«Voo i' Z-ffláE, 

J o 


M 


C n . c , 4 1 


LE 


ÍtP-VJ c 2 


-2 

Ve 2 


r ^ 7 ©^ 

J o 


( 1 ) 


0 -cose) 

d£ = — sin 0 d 9 
2 


y(B) = 2 V, 


(1 + COS0) xp . a 

A.---+> A n smn 0 

0 siné> ^ 


n-í 


With the above, Eq. ( 1 ) becomes 

c mu =- J * A 0 (1 -cos 2 0 ) d 0 - ^ A n (1 - cosQ) sin© sinn 0 d© ( 2 ) 

Note the following definite integráis: 


n -1 


f eos 2 © d© = — 

K 2 

f * sin 2 © d© = — 
J ° 2 


eos© sin 2 © d© = 0 
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f sin9 sin n0 d9 = O 

J o 


for n = 2, 3 


í cos0 sin9 sin 29 d9 = — 

Jo 4 

f cos0 sin0 sin n9 d9 = 0 for n = 3, 4 

J o 


Henee, Eq. (2) becomes: 


o mK =-[nA <} ~^A 0 +^A l -^A 2 ] 


n . A, . 

(Ao + Ai-^) 


4.10 Tile slope of the líft curve is 

ao= °- 65 ~(~ Q39 > = 0.104 per degree 
4 - (-6) 

The slope of the moment coefEcient curve is 

-0.037 - (—0:045) 0 in4 ^ 

mu =---- =8x10 per degree 

4-(-6) 


From Eq. (4.71), 


— + 025 = 


8 x 10" 
0.104 


+ 0.25 
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CHAPTER 5 


5.1 



di x r '= (R di) e 


where e is a unit vector 
perpendicular to the plañe of 
The loop, directed into the page. 


V 



di x r 


T 

\n ^ 0 



r 

4aR 2 


(2?zR)e = 


r -*■ 

-—- e 
2R 


5.2 cU 



Since d i and t are always perpendicular (by inspectíon of the figure). 


_T_ di 

4 n r 2 


dV 


F di x r 


An 


By symmetry, the resultant velocity due to the entire loop must be along the x-axis. Henee, 

T r 2)lR di^ 

I l -r tI I a L 


-| dVcM 0 -[—J, ^-JcpíB» 


1 


Ak (A 2 +R 2 ) 


(27iR) cos9 = 


50 



r r _ r r r 2 

2 (A 2 + R 2 ) 2(A 2 + R 2 ) 372 


5.3 


a = 


1 -i-— (1 + r) 

ttAR 


where ao = 0.1080 per degree = 6.188 per radian 


FromFig. 5.18: 5 = r = 0.054. 


a = 


1 + 


6.188 


(1 + 0.054) 

m 


= 4.91 per rad. 


= 0.0857 per degree 


Cl = a (a - a L =o) = 0.0857 [7 - (-1.3) = |0.712 
C 2 

'-'L 


c„ =- 


?raR 


(1 + 5) = (Q ' 7 1 2) ~ (1.054) = [0.02121 


tt(8) 


5.4 AR= — = -®- = 6.02 
S 170 


At standard sea level, p» = 0.002377 slug/ft 3 


Yoo = 120 mph 


88 fit/sec 


V 60mph J 


= 176 ft/sec 


1 * _ 1 rr\ nnn -• íi — *• z - o n. u¡2 


-P.v- 


(0.002377) (176) = 36.8 Ib/ff 


ao = 0.1033 per degree 
= 5.92 per rad 

L W 2450 


C L 


q ro S q w S (36.8X170) 


0.3916 
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5.92 


a 5 92 

1 + J ^ 1 + 0 1+ -^ (1+CU2) 
ttAR 7t(6.02) 


-4.38 perrad 


= 0.0764 per deg 


C L = a (a - cq=o) 


C L 

a-+ aL=o : 


0.3916 

0.0764 


-3° = 2.121 


5.5 



C 2 


(0-3916) =Q _ 012 g 7 

?r(.64)(6.02) 


D,-=q.S C D . =(36.8X170X0.01267)= 179.3 Ib 


5.6 To be consistent, we wjll use Helmbold’s equations for both the straíght and swept 
wings. 

(a) = 0.1 per degree = 0.1 (57.3) = 5.73 per radian 


ttAR tt(6) 

From Helmbold’s equation for a straigbt wing, Eq. (5.81), 

n- , 

Vl + [a^AR)] 2 +a D /(^AR) 


_5773__ 5.73 

Vl + (Ó.304) 2 + 0.304 ] -349 


4.247 per radian 


(b) From Helmbold’s equation for a swept wing, Eq. (5.82), wbere 
3o eos A = 5.73 eos 45° = 4.05 per radian 


and 


a, cosA 

O 

tzAR 


4.05 

tz(6) 


= 0.215 


we have 
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a 0 eos A 

-Jl + [a o cosA/(7rAR)] 2 +a D cosA/(ttAR) 


_ 4.05 4.05 

Vi+ (0.215) 2 +0.215 1-23785 


3.27 per radian 


Comparing the results of parís (a) and (b), we readily conclude that the effect of wing sweep 
is to reduce the lift si ope. Moreover, the reduction ís substantial. 


5.7 Again, we use Helmbold’s equations. 
(a) ao = 5.73 per radian 
a„ 5.73 


ttAR 7t(3) 


= 0.608 


Vl + [a 0 /(^AR)f+ a 0 /(ttAR) 


5.73 


5.73 


Vi+ (0.608) 2 +0.608 1-778 


(b) a<, eos A = 4.05 


3.222 per radian 


a„ eos A 4.05 


tzAR ;t(3) 


= 0.43 


a„ cosA 


a = 


Vl + [a Q eos A/ (ttAR)] 2 +a o eos A/ (sAR) 


4.05 


4.05 


E.667 


Vi+ (0.43) 2 +0.43 1-5185.- 


ín Problem 5.6, with an aspect ratio of 6, we had 


swept 


3.27 


<W 4 -247 


= 0.77 
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The lift slope for the swepí wing ís only 77% oí that for the straight wing when the aspect 
ratio of both wdngs is 6. 

In Problem 57, with aspect ratio 3, we have 


^swept 
^ straight 


2.667 

3.222 


— 0.83 


Th e lift slope for the swept wing is 83% of that for the straight wing. 

Conclusión : Wing sweep decreases the lift slope. Moreover, wing sweep affects the 
lift slope to a greater degree for higher aspect ratio wíngs than for lower aspect ratio wings. 
This makes some sense, because the lift slope for low aspect ratio wings is already 
considerably reduced just due to the aspect ratio effect. 
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CHAPTER 6 


6.1 V r =—, V 9 = 0, V* = 0 
r" 


Vx V: 


sin# 


c_\ 

da 

c 

~2 


re * 

a 

~dé 


(rsinO) e f 

_d_ 

d<¡> 


1 


r sin 9 


- —{e r (0-0)-r &g 


V d$ J 


4- r sin O e ^ 


o-^íí 

v bq 


i 


r 2 sin 6 


{ 0-0 + 0 } = 


Flow is irrotational. 


6.2 V r =—, V 9 = 0, V*=0 


V V = —— [r 2 (-t)] + — -— (0) +—-— 

r 2 di r' rsin#^# rsin# d<¡> 


V* V = —— + 0 + 0 = 0 + 0 + 0 = 

2 jZ- 

I Oí 


The flow is a physícal possible incompréssible flow. 


6.3 

For the sphere: 


(Cp)-1 



55 


For the cylinder: (C p ) cy i. = 1-4 sin 2 0 

At the top of the sphere: 9 = n/2 s henee 

(Cp) sphere = - 5/4 = 4.25 

For no manometer deflection, (C p ) sp here = (C p ) cy i. 

-1.25 = 1—4 sin 2 0 
sin 2 0 = 0.5625 
sin9 = 0.75 

Henee: 

9^48.6° 

The pressure tap on the cylinder must be located at an angular position 48.6 above or bélow 
the stagnation point. 
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CHAPTER 7 


7.1 p = pRT 

p (7.8)(2116) 


P = 


RT (1716X934) 


= 0.0103 slug/ft j , 


7.2 (a) 



(b) For a calorically perfect gas, c p and Cy are constants, índependent of temperature. 
Henee, we have again 



AJso, at standard sea leve], R = 519°R. Henee 

,6 ñ Ib 


E = 4290 (519) = 2.227 x 10 


slug 


b = 6006 (519) : 


3.117x 10 fc 


ft Ib 


slug 
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-JL-ÍMX287) =10 04. 5 J^ 
r- 1 0.4 kg °K 


r-l 0.4 kg °K 



7.4 


/>. 


RT_ 


4.35 x 10 4 
(287)(245) 


= 0.6186 kg/m 3 


P 

P«> 



vp «7 


P 


fp 


\ l/ r 


Vp w 7 


= 0.6186 


f 3.6 x 10 4 N 

1/1/4 J 

0 5404 ^ 

^4.35 x 10 4 J 


m 


7.5 —; 

Po 


f T \ 


VTJ 


r-i 


T = T 0 


Pl 

kpj 


r-i ( l 

=5or 



0.2857 


E59°fq 


p_ 1.01 X 10 5 = 
RT ~~ (287)(259) 


1.359 kg/nr 


7.6 


pv = RT, henee v = 


RT 

P 
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r ^ 

CN 


\cp. 


RT _ v 

7 ■“ 


Xj-- 


I 


dv 


p 

Note : 1 atm = 1.01 x 10 5 N/m 2 
1 


•c T : 


1 


p (0.2)(1.01 x 10 5 ) 


P.95 x 10' 


-5 m 

N 


Fot an isentropic process: — — 

f/O 

^ 7 


Pl 

V - Pl / 


v.vj 


I.e., piVj Y = p 2 V 2 7 or pv r = constan! = c¡ 




v = 




¿h' 

4» 


= Í (^) , "(p)" 0 "' W =- Í (pv'-) , ' > (p) H ' rt/ ’' = - 1 »p' 

y ' J y 


If&f 


Te - - 


v V.4^ 


2 (--) 

v V ?p7 


7P 


1 


(1.4)(0.2)(1.01 x 10 5 ) 



7.7 


_^ = (L4)Q71Q . 6 o 06 _g_lb_ 

y — 1 (0.4) slug °R 


V 2 V 2 

*ho — h + — -c p T + ——-(6006)(480) + 


(13 00) 2 
2 


3.728 x 10 6 
I slug 


7.8 Let (ho) res - total enthalpy of the reservoir = c p (T 0 ) re s 
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(h 0 ) e = total enthalpy at the exit = Cp T e + 


For an adiabatic flow, ho = constant. Henee 


(bo)res (ho)e 


Cp(T 0 )res Cp T e + 


V e = J 2 c r [fT.>_-T|] = J2(l004.5X1000-600) = {896.4 m/sec 


79 T - P- - (Q-61)(l-01 x 1 0Q 
9 T " p. R" (0.819X287) 


/ ^(r-iVr ( \ / nr 

77 -= — ; T-T.fi = 262.1 ^ 

T„ vp^y vp„y 


= 247.6 °K 


Since the flow ís isentropic, it is also adiabatic. Henee, h 0 - constant 

V 2 V 2 

h«+-f-=h+ — 

2 2 

V- J2(h. -h) + vt = J 2 c,(T_ -T^vl = -J2(1004.5)(262.1 -247.6) + (300)~ 


= 345 m/sec 


V 2 . V 2 
7.10 Pco + p“=p + p-y 


y= f 2 (p^ ~P) ,y 3 . 2 ( 1.01 X 10 3 )(a 61 0-5) + ( - 300 ^2 = 3422 m / S ec 
V K V 0,849 
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7.11 T= T« 


\ t r-TVr 

x 

VPoJ 


; 262.1 


0.3 

0.61. 


: 214 °K 


V = ^2(1004.5X262.1- 214) + (300) 2 = 432 m/sec 


7.12 V = 


_ 2(1.01 x 10^(0.61-03) , 


0.819 


+ (300) =408 m/sec 


432-4087 

% error = I —-— I x 100 = |5.55%| 


432 


7.13 From Eq. (7.53) 


V 2 

h h -= constant 

2 


From Eqs. (7.6b) and (7.9), 


i _ T _ ?RT 

h c p T 

r-i 


From the equatíon of State, 
RT = p/p 

Cómbining Eqs. (1) and (2), 


h= ?- 


y-\ 


V 

\p) 


Henee, Eq. (7.53) can be wiitten as 


r p 
Y~\ \pJ 


V 2 


= constant 


In the limit of y -> °o, Eq. (4) becomes 


0) 

( 2 ) 

(3) 

(4) 
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p. 

p 


V 2 

h -- constant 

2 


or. 


p + Vi p V 2 = constant 

which is Bemoulli’s eqnation. Henee, the energy equation for compressible flow can be 
reduced to Bemoulli’s equation for the case of y —> «>. Henee, the ratio of specific heats for 
incompressible flow is infinite, which of course does not exist ín nature. This is just another 
example of the special inconsistencies associated with the assumption of incompressible 
flow, i.e., constant density flow. which of course does not exist in nature. This is why we 
liave stated earlier in this book that incompressible flow is a mvth . 

As to the qnestion whether Bemoulli’s equation is a statement of Newton’s second 
law or an energy equation, we now see that it is both. For an incompressible flow, the 
application of the fundamental principies of Newton’s second law and the conservation of 
energy are redundant, both leading to the same equation, narnely Bemoulli’s equation. 
However, philosophically this author feels strongly that Bernoulli’s equation is 
fundamentally a statement of Newton’s second law - it is a mechanical equation . This is 
how we derived Bemoulli’s equation in a very straightforward manner ín Chapter 3. For the 
study of inviscid incompressible flow, we need only to apply the fundamental principies of 
mass conservation and Newton’s second law. The principie of conservation of energy is 
redundant and is not needed. 
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CHAPTER 8 









p 0 = 7.824 p = 7.824 (1) = 7.824 ata 
T 0 = 1.8 T =1.8 (230) = 414 °K 
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Since tibe flow is iseníropic, both p u and T 0 are constan.! throughout the flow. Also, in the 
reservoir, M « 0. Henee, the reservoir pressure and temperature are 


p 0 = 7.824 atm 
T 0 = 414 °K 


8.6 From the Standard AMtude 'Pables, at 10,000 ít, 

= 1455.6 lb/ft 2 and T w = 483.04 °R 

T 

From Táble A.1: ForM»=0.82; -^=1.555, —2-®= 1.134 

P. T » 

np 

For M= 1: ^ = 1.893, — =1.2 
' P T 

Since the flow is iseníropic, p 0 = constan! and T 0 = constant. 


„ _ P Po n 

p-Po 

P o P «1 


T T„ 


1.893 


(1.555) (1455.6) = [1196 lb/ft 


1 


T = =— (1.134)(483.04) = [456.5 °R 


T 

8.7 From Table A.2: ^ = 7.72, ^=3.449, —=2.238, 


Pi 


P i 


T, 




9.181, 0-50391 — =0.4601 


Po, 


Henee, 


P 2 


P 2 = — Pi « 7.72 (1) = 7.72 atm) 
Pi 


T, 


T 2 = y-Ti = 2.238 (288) = [644.5 °K 
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(DC1.01 xio-) =1 222 3 

RT, (2B7)(288) 


P 2 = — Pi = 3.449 (1.222) = ^21 kg/ml 
Py 


P^—P. =9.181(1)= [9.181 ,ato| 
Pi 


T 0j = T #i = — Ti = (2.352)(288) = [677.4 °K! (using Table A.1) 


s 2 = Si = * R ^ n ——- ~ (287) i n 0.4601 - 222.8 j ° UlC 


kg D K 


8.8 


Pi_ „ 
Py 


1033. From Table A.2, Uy = 34 -^=2.679, — =12.06 


Pi 


Thus, 


Ti» —T 2 = —í— (1390) = 1518^ 


T, 4 2.679 


From Table A.l, for Mi = 3.0, — = 2.8 


To2 = x ci = _2l t, =2.8 (518.9)= |l453°R 


p 0 , = — p! = (12.06) (1) = [12.06 atm 


Pi 


8 9 = e - ( VV /R = e^ ,99S)/287 = 0.499 

Po, 


From Table A.2: Mi = 2.5 


66 









T 

8.10 From Table A.2: -2- = 2.799 and M 2 = 0.4695 


Henee, 


T 

T 2 = — Ti - 2.799 (480) = 1343,5°R 


a 2 = Va*4)(l716)(1343-5) = 1796.6 ft/sec 


V 2 = M 2 a 2 = (0.4695)(1796.6) = |843.5 ft/sec 


T 

From Table A.1, for M 2 = 0.4695, — = 1.044 

T 


T ^ T 

T 2 * = -4-T 2 = (0.8333X1.044)(1343.5) = 1169°R 

T. T, 


A 2 * = Jtfñz = jmrñWm = 1676 ft/sec 


M 2 *^^5 

‘ a*, 1676 


0.503 


P 1 „ 

'8.11 Is the flow subsonic or supersonic? For soníc ñow, — = --= 1.893, which is 

p 0528 

higher than 1.555. Henee, the flow is subsonic . From Table A.l, for 

. — = 1.555, M = 0.82. 

P 

a = = V(l-4)(287)(288) = 340.2 m/sec 

V = Ma = (0.82)(340.2) = 278.9m/sec 


7712 8 

8.12 The ratio -— = 3.645 is lar ser than 1.893. Henee, the flow is supersonic. This 

2116 * 

means that a normal shock wave exists in front of the nose of the Pitot tube. From Table A.2, 
for 
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P P; _ 7712.8 

p, “ 2116 


3.645,'Mi = 1.56 


a, = = V(1.4)(1716)(519) = 1116.6 ft/sec 


Vi = Mi a 1 = (1.56)(1116.6)= [1742 ft/sec| 


8.13 


p 1.01 x 10 s 
RT “ (287)(288) 


1.22 kg/m 3 


V = | 2 (p °Tp 2= 12(1.55 5 1.0X1.01 x 10 a ) =3Q3lIl/sec incORRECT 
\ p V 1-22 


% error = 


(b) p = 


303-278.9 _ 
278.9 

p 2116 


!8.69°/d 


RT (1716)(519) 


0.002376 slug/ft 3 


y = 1 2 ^°.. . ? 2 = =2170.5 ft/sec TNCORRECT 

. i p V 0.002376 


n/ 2170.5-1742 _ F ^ 71 

% error =---124.6%| 


1742 


8.14 


p2_ n 2 r jY r + i+2y M;-2 r _ i r + 2/ Mr 

p, x+i 1 x + i x+i 


(i) 


Po, 
P 2 



i+[(r-i)/2]M, 

2 >MÍ-(r-i)/2) 


( 2 ) 


Working vjitii the expression inside the parenthesis of Eq. (2): 
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i + ¿r±M’.i + ^± 


1 + - 


r — 1 


Mr 


y Mí ~(y-l )/2 


1 + ( Y -1) 


1 + 1 ——-1 M? 


27 Mí-(r-i) 


= 1 + (y - 1 ) 


2+p-i) Mí 
4 y Mí-20 — 1) 


4 y MÍ-2p-l) + 2p-l) + p-l) 2 M? 
4y Mí -2 0-1) 


4y Mí + p 2 ~2y + l) Mí _ O 2 +27 + 1 ) MÍ 
4/ MÍ -2 0-1) 4 y Mí -2 p-1) 

0 +l) 2 MÍ 
4 y Mí-2 0-1) 

Combining Eqs. (4), (2), and (1), we have: 


Po, Po ; . p 2 _ 

0 + l) 3 MÍ j 

7 

7-1 

I -7 + 27 MÍ 

Pl p 2 Pl 

_4y Mí-2 0-1)„ 


7+1 


which is Eq. (8.80) 


8.15 At 80,000 ft., T ffi = 389.99°R 

V ao = 2112 = 3097.6 ft/sec 

a, = (1.4)(1716)(3 89.99) = 967.9 ft/sec 


M*,= 


3097.6 

967.9 


= 3.2 


From Appendix A: 

ForMco = 3.2, =3.048 

T 

fio 

T 0 = 3.048 T» = 3.048 (389.99)= 1188.7 °R 
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Since 0°F = 460°R, the 


T 0 = 728.7°F 


8.16 



1.13 

0.1 


11.3. 


From Appendix B, M® = 




8.17 The temperatme at the stagnation point is the total temperarme in the freestream, 
because the total temperature is constant across the normal shock. From Eq. (8.40), 

. — = l + = 1 + -^—^ (36) 2 = 260.2 

2 2 


Since Tco = 300 K, we have 


T 0 - (260.2X300)- [78 ? 060K 


This is an ungodly high temperature. It is also íncoxrect, because long before the aii would 
reach this temperature, ít would chemically dis-sociate and ionize. In such a chemically 
reacting gas, the specific heats are not constant which means that Eq. (8,40) is not valid for 
such a chemically reacting flow. In realíty, the temperature at the stagnation point on the 
Apollo was cióse to 11,000 K, much lower than our estímate above, but still plenty high. Air 
at 11,000 K is a partially ionízed plasma. For the analysis of high temperature, chemically 
reacting flows, techníques much different than those discussed in this book must be used. 
See for example Anderson, Modem Compressíble Flow , 2 nd ed. ? McGraw-Hill, 1990, or 
Anderson, Hypersonic and High Temperature Gas Dynamics , McGraw-Hill, 1989, reprinted 
by the American Institute of Aeronautics and Astronautics, 2000. 


8.18 Use Eq. (8.40) 


A = i + 2^J-m; 

T 2 


For T 0 = 11,000 K, T® = 300 K, and M® = 36, this equation becomes: 


l 1 , Q 00 = 1+ y__l 36)2 
300 2 


35.67 = 648 y - 648 
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or, 


683.7 

648 


1.055 


In order to use Eq. (8.40) to estímate a reasonably valid stagnation temperature for the 
Apollo, we have to use an “effective gamma” of 1.055. To double check this, retum to Eq. 
(8.40), insert y = 1.055, and calcúlate T 0 . 


A =]+ üz1 mí 

T 2 


1 + L ° 55 1 (36) 2 = 36.64 


or. 


T„ = 36.64 = 36.64 (300)= [11,000 K 
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CHAPTER 9 


9.1 



h ® 559 Tan p = 559 Tan 41.8° 


fa = 500ffl 


9.2 M B] = Mi sin (3 = (4.0) sin 30° = 2 


From Table A.2, for M #i =2: — = 4 . 5 ; — 1.687, — ; — 0.7209, M nj 0.5774 


Pi 


Po, 


P 2 = pi = ( 4 . 5 ) (2.65 x 10 4 ) = 11.193 x 10" N/mj 
Pi 


T 2 = — T] = (1.687)(223.3)= l376.7°K 


From the 0-p-M diagram: 0 = 17.7° 


M 2= ^_ = —<*™-=£ñ 


sin(/? — ff) sin(30 — 17.7) 

n T 

From Table A.1, for Mi = 4: = 151.8, = 4.2 

Pi T i 


xi /_ A 


p oj = pi = (0.7209)(151.8)(2.65 x 10 4 ) = l2.9xlQ°Wnf 


Po, Pt 


T 0j = T 0 _ = -2L Ti - (4.2)(223.3) = [937.9°K 
Ti 
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S 2 — si = -R £ n 


Po, 


- (287) 0.7209 : 


93.9 


joule 

WK 


9.3 Consider an oblique shock.. For such a case, 



Depends on actual Mach Depends on notmal 

number behind the shock Mach number upstream 

M?, not M„ of the shock, M„ . 


In the derívation of Eq. (8.80), we related M 2 directly to Mi thxough Eq. (8.78). This holds 
only for a norm al shock. If we wish to use Eq. (8.78) for an oblique shock, then both M 2 
and M] in Eq. (8.78) are replaced by M n; and M . However, in Eq. (1) above, p 0j /p 2 

Depends on M 2 , not M n;¡ . Because Eq. (8.78) does not relate M 2 to Mj for an oblique shock 
(it relates M n to M ), then Eq. (8.78) cannot be used for the derívation of p 0j / p, for an 

oblique shock. Therefore, the derívation of Eq. (8.80) holds only for a normal shock. It can 
not be used for an oblique shock, even with Mj replaced by . On the other hand, 

D T 

S 2 -Si=c p ^n — - R f n — 

Pi T i 


wherep?/pi and T-/T 1 for an oblique shock depend onlv on M n . Since = e^Y*^then 

Po, 

P 

clearly —— depends only on M n . For these reasons. when using Table A.2 to determine 
Po, 

chanees across an oblique shock, using M„ , the total pressure ratio j s a valid column, 

Po, 

but the column giving ■^ L * is not valid. 

Pi 


9.4 To CORRECTLY calcúlate Poj : 

M ai = Mi sin (3 = 3 sin 36.87° = 1.8 
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From Table A.2, for M = 1.8: ^ = 0.8127 

P„ 

r'o» 


From Table A.l, for Mi = 3: 


= 36.73 


p| = (0.8127)(36.73)(1) = 129.85 atm 


P 0l Pi 

(b) The INCORRECT calcuiation of p 0j would be as follows: 

Po, 


From Table A.2, for M n =1.8: 

7 "i 


; 4.67 


p —— pi = 4.67 (1 atm) = 4.67 atm. Totally WRONG 
° 2 Pi 


% eiror = x 100 = 539% - a terríbly large error. 

4.67 


9.5 



From tibe 0-p-M diagrama |P = 46° 

M ni =Mi sin p = 2.5. sin 46° = 1.8 


j 1 

From Table A.2, for M„ = 1.8, =3.613, -f = 1.532, M Jlj =0.6165 

' Pi h 


p 2 = — pi = 3.613 (1 atm) = |3.613 atm| 
Pi 
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T] = (1.532)(30Q) ~ 459.6°K 


M 2 = 


M, 

n 2 

Sin(/?-<9) 


0.6165 

Sin(46-22.5) 


1.546 


9.6 From the 9-p-M diagram, shock detachment occurs when a > 28.7°. When a = 9 = 
28.7°, (} = 64.5°. 


M n) = Mi sin p = 2.4 sin 64.5° = 2.17 


From Table A.2, for M n = 2.17: = 5.327 

7 Dj 

P» 


Pi 


Pmax = ~ Pi “ 5.327 (1 atoa) = p.327 atmj 
P. ~ ~ 


and the máximum pressure occurs when a — |28.7° 



From the 9-P-M diagram: P = 48° 

M ni = M} sin P = 3.5 sin 48° = 2.60 

From Table A.2: -^*-=0.4601. =0.5039, 

n 2 J 
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m 2 


M. 


0.5039 


Smtf-ff) Sin(48 - 30.2) 


1.648 


From Table A.2', for M 2 = 1.648; = 0.876 

Po, 


From Table A.1, for M = 3.5: ^=- = 76.27 

P» 


Po, 


Po, Po, Po 


Po, Po. P« 


p»= (0.876)(0.4601)(76.27)(0.5) = 115.37 atm 


9.8 From Table A.1, for Mi = 4, —= 151.8 

P. 


Henee, p 0> 


p, = 151.8(1 atm) = 151.8 atm. 
P. 


a) Mi = 4 


From Table A.2, for M] =4: ^--0.1388 

Po, 


Po, 


Po 2 


Po 2 


Í22- p = 0.1388 (151.8) = 21.07 atm 

Po, 


Loss in total pressnre = p 0j - p 0j 


151.8 —21.07 = 


130.7 atrd 


b) 
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From the 9-(3-M diagram, 

P = 38.7° 

M n _ = M, sin p = 4 sin 38.7° =2.5 


From Table A.2, for M n = 2.5: ^-= 0.499, M =0.513 

"1 „ n 2 


M 2 = 


M 


0.513 


$in(/? — <?) sin(3 8.7-25.3) 


= 2.21 


From Table A.2, for M 2 = 2.21: ^ = 0.6236 

Po, 


Po = — — — Pi = (0.6236X0.499X151.8)(1 atm) = 47.24 atm 
J Po, Po, P] 


Loss in total pressure = p Di - p 0i = 151.8 - 47.24 = ¡1Ó4.6 atro 

c) 



From parí (b) above, Mo ' 2.21, -2— = 0.499. 

' ' Po, 


From tibe p-0-M diagram: p 2 = 47.3 c 
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For the second shock: = M 2 sin P 2 = 2.21 sin 47.3° - 1.624 


From Table A.2, for M„ = 1.624: ^ = 0.8877, M nj = 0.6625 

■ . 2 Po, 


M 3 


M. 


0.6625 


sin(/? 2 —0 2 ) sin(47,3—20) 


P 0J 


- 1.444 


From Table A.2. for M3 = 1.444: - 0.947 

Po, 


p Eli. E^ p¡ = (0.947)(0.8877)(0.499)(151.8) 

" Po, Po ; Po, P. 


p„ = 63.68 atm 

J- 


Loss ín total pressnre — p o - p 0 = 151.8 - 63,68 - [88.1 atm] 


CONCLUSION : To decrease a supersonic flow to subscmic speeds via a shock svstem, a 
series of oblique shocks followed by a normal shock yields a smaller total pressure loss than 
a normal shock by itself. Henee, a svstem of oblique shocks, followed by a normal shock is a 
more efficíent means of slowing a supersonic flow to subsonic speeds than a single normal 
shock itself. 


9.9 
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From the 9-P-M diagram, p l = 34.2' 
M n| = Mi sin p¡ 


= (3.2) sin 34.2° = 1.8 

From Table A.2; for M n = 1.8: —=3.613, — = 1.532. 

P. T, 


M„ =0.6165 


M 2 = 




0.6165 


sin(j5, - <?,) sin(34.2 — 8.2) 


: 2.24 


For the Reflected Shock : 

From the 0-P-M diagram, for M¿ = 2.24 and 9=18.2°: p 2 = 44° 
M„ 2 = M 2 sin p 2 = 2.24 sin 44° = 1.56 


HT 

From Table A.2. for M n = 1.56: ^ = 2 .673, = 1.361, 

- n ? 7 ’T' 7 

p 2 t 2 


=0.6809 


M 3 = 


M. 


0.6809 


sin(/? 2 - ff) sin(44 -18.2) 


1-56| Note : The fact that M 3 and M nj are 
equal is just a coincideñce. 


Q = p 2 -0 = 44-18.2 = ¡25.8° 


p 3 = ill ib p] = (2.673)(3.613)(1 atm) = ¡9.66 atm 
P 2 Pi " " 

j 3 = 5.2* Ti = (1.361X1.532X520) = jl084°R 

t¡ 
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9.10 



From Table A.3: ForMi =2. v, =26.38° 

v 2 = 0 + Y, == 23.38° + 26.38° = 49.76° 


Henee, 




T 

From Table A.l,forM,= 2: ^ = 7.824, 


Pi 


T, 


1.8 


For M 2 = 3: ^- = 36.73, -£- = 2.8 


However: p„ = p„ and T = T . Thus 

r °1 r O) Oj O; 


p 2 = p, = (-J—| (7.824X0.7) = 10-149 atm| 

P 0 , P, v 36.73. 


T 2 = Tu ^ T, - ÍTI (1.8X630) = twsg 

i, 0 ii V2.ó, 


p, - «Uggllgl . ^ x 1¡r 


RT, (1716X405) 


Pc 2 = Po, = — Pi = (7.824X0.7) = 15.477 atm¡ 
P. 
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T 


T_ = X. = — T¡ = (1.8)(630) = [l 134°R| 


From Table A.3: for = 2, pj = 30° 
For M2 = 3, p 2 = 19.47 
Referenced to the upstreara direction: 


Angle of forward Mach line = pi = 30 


Angle of rearward Mach line = u 2 -9= 19.47 - 2338° = [-3.91 


Note: The rearward Mach line is below the upstream. direction for tbis problem. 


9.11 From Table A.1 T for M, = 1.58: ^ =4.127 

Pi 


= Üí». = ÍL = (4.127) í—1—1 =31.6 
p 2 p 2 pj p 2 V 0.13067 


From Table A.1. for ^ = 31.6, M 2 = 2.9 
P2 


From Table A.3, for Mi = 1.58; vj = 14.27, for M 2 - 2.9: v^- 47.79 


9 = v 2 -v, = 47.79- 14.27 = 


B3.524 


9.12 
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From the O-p-M diagram: 


For M, = 3 and 0 = 30.6°, p = 53.1 ° 
M n| — sin p = 3 sin 53.1 = 2.4 


Pz _ C £ X ^ j\ a Po, 


From Table A.2, for M = 2.4: = 6.553, ^ = 2.04, ^ = 0.541, M n = 0.531 

Pl X Po, 


M n 

M 2 =- 5 


05231 


sin(y9 - 6) sin(53.1 - 30.6) 

From Table A.3: For M 2 = 1.37, v 2 = 8.128 
v 5 = 8.128+ 30.6 = 38.73° 


1.37 


From Table A.3: For v 3 = 38.73°, M 3 = |2AS 

P. - X 


From Table A. 1: For Mi = 3, 


36.73. =2.8 


Pi 


T 

For M 3 = 2.48,: ^=16.56, ^=2.23 

P 3 T 3 


P3 : 


Ps_Po^ Po 
Po, Po 2 Po 


Pl 


■Pl 


1656. 


(1)(0.5401)(36.73)(1 atm)= [120 atml 



X 

X 





(1)(1)(2.8)(285) = |357.8°K 


Clearly, p 3 p 1? T 3 ^ Ti, and M 3 + M. Why? Because tbere is an entropy increase across 
shock wave, which permanently alters the thermodynamic State of the original flow, even 
after it is brought back to its original direction. 
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(a) For Mi = 2.6 and 0 = 5°, p = 26.5° 

M n| =Mi sin p = 2.6 sin 26.5° = 1.16 


From Table A.2: —=1.403 
Pi 


From Table A.l, for Mi = 2.6: — = 19.95 

P, 


From Table A.3, for Mi = 2.6: Vi = 41.41 ° 


v 2 = vj+ 0 = 41.41 +5° = 46.41°-+M 2 = 2.83 


From Table A.l, for M 2 = 2.83: — = 28.4 

P 2 


Pl = Ja = (0.0352)(1)(19.95) = 0.7022 

Pl Po, Po, Pi 



(P ; -P 2 ) c c «s « 
q„ c (1) 


(P 3 ~P Z ) 

q« 


eos a 


q<» ” qj ~ Vi p) Vi 2 = 


r PiA v f 

2 r p. 


Y P.1 V T 2 _ Y P,M| 
2 af 2 
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_ 2 (p 3 ~P2) 

r Pi M 2 


eos a 


X M? 


P 3 P 2 
VPi P,/ 


eos a 


(1.4)(2.6 ) 2 


(1.403 - 0.7022) eos 5 o = j0.148| 


c d = 


2 í P 3 Vi 


7 M 7 Vp 2 p,7 


sin a _ . sin 5° rX. 

sin a = c e -= 0.148 —— = (0.01291 


cosa 


cos5° 


(b) For Mi = 2.6 and 9 = 1 5°, P = 35.9° 

M„ i = M, sin (3 = 2.6 sin 35.9° = 1.525 

From Table A.2: — = 2.529 

P. 


From Table A l, for M, = 2.6: ^ = 19.95 

Pi 

From Table A.3, for Mj - 2.6: v, = 41.41° 

v 2 = V! + 0 = 41.41 + 15 = 56.41°M 2 = 3.37 

From Table A.l. for M 2 = 3.37: = 63.33 

P 2 


P 2 _ P 2 P», 

Pl Po, Po, Pt 


*-(-U 

p, V 63.337 


(1X19.95) = 0.315 


£ Y M? 


\ 


Pl_P2 
VP, Pi 7 


eos a = 


(1.4)(2.6 ) 2 


(2.529 - 0.315) eos 15° = |0.452| 


sin a sin 15° 

c d = c, -=0.452 •—-— =10.121 


cosa 


cosí5° 


(c) 


For Mi = 2.6 and 0 = 30°, p = 59.3° 
M n - Mj sin p = 2.6 sin 59.3° = 2.24 
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= 5.687, ^ = 19.95, v, = 41.41° 
Pi P. 


v 2 = V! + 0 = 41.41 + 30 = 71.41 0 -> M 2 = 4.46 


^ = 275.25 
P2 


= Esl 

Pl Po 2 Po, Pl 


í-í- 

V 275.25. 


(1)(19.95) = 0.0725 


Y 


íp> 

vp, 


p/ 

vJ 


eos a = 


(l-4)(2.6) 


(5.687- 


0.0725) = |1.19| 


c„*U9 


sin 3 0‘ 



For región 2: 


v, = 49.76° 


v 2 = vj + 0 = 49.76° + 5° = 54.76° M 2 = 3.27 

For Mi = 3: ^ =36.73: 

Pi 


For M 2 = 3.27, = 54.76 

P 2 

For región 3 : 

v 3 = v 2 + 0 = 54.76° + 20° = 74.76° M 3 = 4.78 

ForM 3 = 4.78: ^ =407.83 
P3 


For región 4 : 


Mi = 3 and 9 = 25° -» p = 44° 


= Mj sin p = 3 sin 44 = 2.08 


^ =.4.881, M n = 0.5643, and = 0.6835 
Pi 4 Po, 


M 4 = 



$in(/?— ff) 


0.5643 
sin(44 - 25) 


1.733. 


Thus, 


v 5 = 18.69, ^-=5.165 

P4 

For región 5 : 

v 5 - v 4 + 9 = 18.69° + 20° = 38.69° -> M 5 = 2.48 
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^ = 16.56 
Ps 


Pressure ratios 


P 2 _ P 2 Pq, Pq, _ 
Pl Po 2 Po, Pl 


- 2 —] (1X36.73) = 0.6707 
54.76/ 


Pi = P^Pl = P2.l2_Z£a.Z^ = (0.6707) í ——) (1X54.76) = 0.09 
Pi Pi P 2 P. P 0; P 0i P 2 V 407.837 


=4.881 

Pi 


P 5 _ Ps Po; Ppj Po, 
Pl Po, Po 4 Po, Pl 


-¿-](1)(0.6S35)(36.73)= 1.516 
16.5 o/ 


Let £ = length of each face of the diamond wedge. 

L-' = p4 l eos 25° + p5 £ eos 5 o - P2 í eos 5 o - p 3 eos 25° 
L' = (p 4 -p 3 ) l eos 25° + (p 5 -P 2 ) i eos 5 o 




2 t 
j M 2 c 


r/ 



r 

A 

“ 


cos25 D -b 

p 5 

_P 2 

cos5° 

LVp, 

Pi7 


Ipi 

Pj 



——- - [(4.881-0.09) eos 25°+ (1.516-0.6707) eos 5 o ] 

(1.4X3 ) 2 c 


9 = 0.823 -. 

c 


However, 

. c /2 
i 


= eos 10° - = ■■■■■ 1 -- - = 0.5077 

c 2 eos 10 ° 


9 = (0.823)(0.5077) = 10.4181 
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D' = p 4 í sin 25° + ps £ sin 5 o - p 2 £ sin 5° - p 3 £ sin 25° 
D' = (p 4 — P 3 ) £ sin 25° + (ps - P 2 ) £ sin 5 o 


Cd = 


D’ 

q«s 


D’ 



2 í 
y M 2 c 


rr ^ 

r 

A 1 

£i_pi 

| sin 25°+ 

.Pj. 

-—] sin5° 

.Vpl PiJ 


^P) 

P J 


2 £ 
(l-4)(3 ) 2 c 


[(4.881 - 0.09) sin 25° + (1.516 - 0.6707) sin 5 o ] 


c d = 0.333 - = 0.333 (0.5077)= |0.169| 
c 


9.15 The máximum expansión would correspond to M 2 —> °o. From Eq. (9.42) in the text, 


lim v 2 = lim 


tan- üE. 

\y-l Vr+11 


(M 2 - 1 ) - tan " 1 ^M 2 -1 


T'v / Í2 ^ ÍVÍ 2 —^ 


= L±1 £ .£ . ( 2£tl - i) — =2.277 rad - 130.45° 

\y-\ 2 2 lr-1 2 


Since, for Mi = 1, vi = 0, then 


0 = V 2 - vi = 130.45 - 0 = 1130.45' 
max 



/30.4S 


O 



9.16 For the cylinder, with c d based on frontal area. 


(D')cvi = q w Scd = q m d(l )/(4/3) = -j (d) q„ 

For the dimensional wedge airfoil, refemng to Figure 9.27. 
(D')w = (P 2 - P3> t 


Henee, 


(D ')* _ I (d) q ” 

(D% (p 2 -p 3 )t 

y 2 

However, t = d and q„ = —■ pi Mi 
Thus, 


(D%, 

(D*)w 


3V2y 



I r M? 

J _ 

P2 P3 

vpi p,y 


To calcúlate p 2 /pi, we have, for Mi = 5 and 0 = 5 o , p = 15.1°. 
M M =M! sin p = 5 sin (15.1°) = 1.303 

From Appendix B, for Mnj — 1.302, - 1.805. Also, 

P. 


M 2 = 


M 


n,2 


sin(/? — 0) 


0.786 

sin(l5J - 5) 


= 4.48. 


To calcúlate —, the flow is expanded through an angle of 1 0 o . From Table C, fox M 2 
P. 

4.48, v 2 = 71.83 (nearest entry). 

v 3 = v 2 +9 = 71.83+ 10 = 81.38° 


Henee, M 3 = 5.6 (nearest entry) 
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From Appendix A: For Mi = 5, —= 529.1 

Pj 


ForM 3 = 5.6, ^ =1037 
P 3 


p 

From Appendix B: For M n = 1.303, =0.9794 

' Po, 

Thus, 


- _ P? p°, P°z P»L = f_L_ ) (1X0.9794)(529.1) = 0.5 
P, P„, P„, Po, P, V1037/ V 

Henee, 

(D’)„ \r m; f (1.4)(5) ! 

( D ). fp 1 _p,'| (1.805-0.5) 

Vpj P J 

Note : This is why we try to avoid blunt leading edges on supersonic vehicles. (However, at 
hypersonic speeds, blunt leading edges are necessary to reduce the aerodynamic heating.) 


17.a 


The flow direction downstream of the leading edge is given by line ab. The flow direction is 
below the horizontal (below the direction of M») because lift is produced on the fíat píate, 
and due to overall momentum consideratíons, the downstream flow must be inclined slightly 
downward. Also, line ab is a slip line; the entropy in región 4 is different than in región 5 
because the flows over the top and bottom of the píate have gone through shock waves of 
different strengths. The boundaiy condition that must hold across the slip line is constant 
pressure, i.e., p 4 = p 5 . It is this boundary condition that fixes the strengths of the expansión 
wave and the shock wave at the trailing edge. 

To calcúlate the trailing edge shock and expansión waves, and the flow direction 
downstream, use the followíng íterative approach: 

1. Assume a valué for <j>. 

2. Calcúlate the strength of the trailing edge shock for the local deflection angle (a- 
<j>).- This gives, among other quantities, a valué of p 4 . 

3. Calcúlate the strength of the trailing edge expansión wave for a local expansión 
angle of (a.-<j>). This gives a valué for ps- 

4. Compare p 4 and p 5 ffom steps 3 and 4. If they are different, assume a new valué 
of<j>. 

5. Repeat steps 24 until p 4 = p 5 . When this condition is satisfíed, the iteration has 
converged, and the trailing edge flow is now determined. 


CHAPTER 10 


10.1 From Table A. 1 , for AJA* = 2.193, 


^ = 12.5, i- = 2.058. 
P< T c 


For isentropic flow, T 0 = constant and p 0 = constant. Henee, 


Po e =Po = [5 atm| , and = T 0 = 




Te= -T 0 = í_i— ] (520)= \252J°l 
T V2.0582 e -- 


p e =^- = (0-4X2116) - |o oo 195 slug/f 
RT e (1716X252.7) 1 - *~ 


I (1.4)(1716)(252.7) = 779.2 ft/sec 


u e = M e a e = (2.3)(779.2) ■ 


10.2 — =-— ”3.182. From Table A.l,weseethat 

p e 0.3143 


and AJA* 


103 Ahead of the normal shock in front of the Pitot tnbe. 


p 0i = p 0 = 2.02 x 10 5 N/m 2 


P^ = 852X101 = 0.4416 
p„ 2.02 x 10 5 

r °l 
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From Table A.2: M e — 2.65 


From Table A. 1: A e /A* = 13.036 


1A . ' * * A *. n Pn _ (5)(2116) 

10.4 m p n ,p 0 (1716)(520) 


186 


slug 

ñ 3 


p* = £l p o = (0.634)(0.01186) = 0.007519 slug/ft 3 
Po 


T* = _L_ To = (0.833)(520) = 433.2 D R 
T 


u* = a* = -704X1716X4332) = 1020 ft/sec 
m = p*u*A* = (0.007519X1020) 


U-U 

0.213 SlUE 

VI447 

sec 


10.5 m = p*u*A* 

u* = 0T* and p* =■ i. 


Henee, 



A* V^ T * 


p*A* 

RT* 



Since, M* = 1, tben 




= l + 'LJl M * 2 = 
2 


y +1 

2 _ 


P a 

P* 



Thus, 
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10.6 p 0 = 5 ato = 5(2116)= 10580 lb/fi 2 
A* = 4/144 = 0.02778 fe 


¿ = (105809X002778) (1.4) f 2 V _ slug 

V520 y (1716) V 24J sec | 

which is the same as obtained in 
Problem 10.4 


10.7 
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First, check to see íf the flow is soníc atthe throat. 


=1.056 

p c 0.947 


From Table A., for — = 1.056: M e = 0.28 and A e /A* = 2.166 
Pe 


Since A* = 1 616 < =2.166, tben A, > A*. The throat size is larger than that for sonic 

A, A* 

flow, henee the throat Mach nnmber ? M t , is snbsoníc. 

Aí_ = Aí, As. = _JL_ (2.166) =1.34 
A* A e A* 1.616 



1 0.8 Note : The equation for m given in Problem 10 5 can not be used here because the 
flow is not chokedj Le. ? the throat Mach number is not sonic. 


» 

m = Pe Ae U e 

From Table A.l.for — —=1.056: M e = 0.28, —f- =1.016 
Pe J -= 

Te = To/1.016 = 288/1.016 = 283.5°K 

p - (°-947)(1.01 x 10 5 ) = , 1?6 kg/m 3 
Pt RT e (287)(2835) 

A = JW' = a/( 1.4)(287)(2835) = 337.5 m/sec 

Ue = Me a e = (0.28)(337.5) = 94.5 m/sec 
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A e = A t 



VAJ 


(0.3)(1.616) = 0.4848 m 2 


m = p c A e u e = (1.176)(0.4848)(94.5) = |53.88 kg/se C 


10.9 (a) ~ = —i— = 1.064. 

p e 0.94 


From Table A. 1: M e = 0.3 and AJA* = 2.035. —- 

A* A. A* 


(2.035) = 1.33. 


A ‘ A - A = (JL] 

153J 

Since A t > A*, then the flow is completely subsonic. No shock wave exists. Henee, from 
P» 

P e 


Table A. 1, for = 1.064, M« = 0.3| 



1 

0.886 


1.129. 


From Table A.1, for ^ = 1.129: M e = 0.42 and = 1.539. 


A_ = A¡_ A. 

A* A e A* 


—^—1 (1.529) = 0.999 « 1.0. 
.1.537 v 


Henee, A t = A*, and theflowisprecisely sonic at the throat. It ís subsonic everywhere else. 


Henee, from the above M e = 0.42 


(c) From tlie ábove results, clearly when p e is reduced below 0.866 atm, sonic flow 
will occur at the throat, and the nozzle will be choked. Since p e = 0.75 atm is far above the 
supersonic exit pressure, we suspect that a normal shock wave exists within the nozzle. Note 
that, if we run the same calculation as in parts (a) and (b) above, we find: 



From Table A.1, for — - 1.333, we have 
P e 


—s- = 1.127 
A* 
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A t _ Al Al. = (—— ] (1 127) = 0.7366. Since it is impossible for A t < A*, then 
A* A e A* U537 

clearly the flow can not be completely isentropic. There mnst be a shock wave mside the 
nozzle, wíth a consequent change in both p 0 and A* across the shock. Henee, e a ove 
calculation is meaningless. bastead, set up the following tiial-and-error process as follows: 

Assnme a normal shock exists insíde the nozzle, say at a locatíon where A 2 /A t - 
1.024. Let: 


Ai* = sonic throat area for the flow ahead of the shock. 
^ 2 * — sonic throat area for the flow behind the shock. 
p o = total pressure for the flow ahead of shock. 

p o = total pressure for the flow behind shock. 


Shock 



Note that p 02 < p 0( A 2 *>Ai* 

which comes from the shock wave theory discussed in the text 
Kev equation : 
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( 1 ) 


p= = ~ —p OI 

Po, Po, 

To find the valúes of the ratios m Eq. (1): 

From Table A.l for A 2 /Aj * = 1.204: M] = 1.54 

From Table A.2 for M, = 1.54: M 2 = 0.6874, ^ = 0.9166 

Po. 


From Table A.1, for M 2 = 0.6874: —= 1.1018 

A * 


~~~ = A. Al A_ = (l.53)í —i— j (1.1018) = 1.4 

A 2 * A t A 2 A 2 * U.2047 


From Table A.l. for = 1.4: M e = 0.47, = 1.163 

A, * o 


Retuming to Eq. (1): 
P 


P e 


.(-i: 

P„, P “ UJ6Í 


(0.9166)(1 atm) = 0.788 atm. 


Ibis is slíghtly higher than the given p e = 0.75. Henee, move the shock wave slightly 
downstream. 

Assume A 2 /A t = 1.301 

From Table A.l: Mi = 1.66 


From Table A.l, for Mi = 1.66: —— = 0.872, M 2 = 0.6512 

Po, 


From Table A.l, for M 2 = 0.6512: —— = 1.1356 

A 2 * • 


A A A A i 

—^ = ^ =(1.53) T ~ T |(1.1356) =1.335 


A 2 * A t A 2 A 2 


1.301 
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From Table A.l, for= 1.335: M e = 0.50, ^- = 1.1862 
A * p. 


From Eq. (1): 

_ Pe Po 


Pe = 


Po, Po, 

A. 


Po, 


( 1 


V 1.1862 


(0.872)(1 atm) = 0.735 atm. 


q r 

Interpólate: 4Ü = 1.301 — (1.301 — 1.204)---- = 1.274 

v A, 0.788-0.735 


Thus, Assume A 2 IA l = 1.274 

From Table A.l: Mj = 1.63 


From Table A.2: M 2 = 0.6596, ^ = 0.8838 

Po, . 

From Table A.l: - 2 . = 1.1265 
A 7 * 


A e = A e =( 1 .53)í_l_| (1.1265) = 1.353 


a 2 * a, a 2 a 2 


1.274. 


From Table A.1: M e = 0.49, — =1.178 

Pe 


p c = —— p o =1 —í—] (0.8838)(1 atm) = 0.75 atm 


_ Pe Po 


1 


Po, Po, ’ ''1378/ 

Henee. p e calculated agrees with p e given. Thus, 
M e = ÓÁ9| 


(d) 



1 atm 
0.154 atm 


6.49-' 


From Table A.l: —- = 1.53. which ís precisely the given area ratio of the nozzle. Henee. 
A* 


for this case, we have a completely isentropic expansión, where. 
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M e ~ 1-88 


10.10 From the 0-p-M diagram, for 0 = 20° and p = 41.8°, we have M) = 2.6. From Table 
A.l, 


A 5 




10.11 From Table A. 1. for ^ = 6.79, M c = 3.5 

' A* 


From Table A.2, for M e = 3.5: —— = 0.2129 

Po, 



( 1 ) 


-- (1.448) = 

v 0.2129/ 1 

6.8 atm 


10.12 From Table AT, for M e = 2.8: ^=27.14,-^=2.568 

P* T e 


At standard sea level: p = 2116 lb/ft 2 , T = 519°R 



Po = — =--= 0.251 slug/ft 3 

RT 0 (1716)(1333) 

p* = (0.6339)(0.0251) = 0.0159 slug/ft 3 - 
T* = 0.833 (1333) = 1110°R 

a* = V;RT* = ^J(lÁy(í716}^ñO) = 1633 ft/sec = u* 
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m = p*u*A] 


A, 


* -~ 


m 


1 


p*-u* (0.0159X1633) 

From Table A-l: AJA* = 3.5 


= k).0385 ft 


A c = A* = (3.5X0.0385) = |Q.1348~1 4 


From Eq. (10.38) in text: 


A t; _ A 2 * _ P 0¡ 
A,, A, * p. 


From Table A.2: for M c — 2.8i — 0.3895 

Po, 


A = A, 

h m 


.Po J 


= af 


'p3 


VPoJ 


(0.0385)í—— 
v V 0.3 895. 


0.0988 ft 


10.13 m = p*u*A* 

Also,R^M = ^ = 519.6^ 

16 kg K 


( 1 ) 


P* = — Po = 


T* * 

T* = T = 
T„ 




y + V 


( 2 ^ 


2V 


y + 1. 


RT e v 2.2/ (519.6)(3600) 

(3600) = 3273 K 


= 3.319 xl0’ 7 p o 


u* = a* = V?Rf* = 40 -2X519.6X3273) = 1428.6 m/sec 

kg 


Henee, firom Eq. (1), with m - 287.2 


sec 


287.2 = (3319 x 10‘ 7 p o )(1428.6)(0.2) 


or. 
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Po = 


287.2 


(3.319x10~ 7 )(1428.6)(0.2) 


= 3.029 x 10 6 


N 


m - 


or. 


3.029 x 10 6 
1.01 x 10 5 


= 30 atid 


10.14 We assutne the flow velocity is low at the diffuser exit; henee the total pressure at the 

exit is 1 atm. From Appendix B, for M = 3. 0.3283. 

Po, 


no = JbAi- = 1.2 

Po 2 /P 0l 


Pb _ 1 O P°2 _ ■ 


^ = 1.2 = 1.2 (0.3283) = 0.394 

Po Po. 


Po ! 


Pb , l _ 
0394 0.394 


E.54 atad 
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CHAPTER11 


11.1 u = Ét = y. + (7Q) - eos (2 tix) 

^ Jl-Mi 


d<j> __ 70 




(2^1 tí) 


^-í.-ryt-M^y 


sin (2-nx) 


- - 140 7 t e 2,c 7'-M-y gjjj (2 tdí) 


a«= = VC^X 1716 )^ 19 ) = 1116 - 6 ñ/sec 


M*, = — = = 0.6269 

a_ 1116.6 


HiuS; at (x,y) = (0.2. 0.2) 

u= 700 + e -2yo.779)(02> CQS p 7 i(. 2 )] = 765.6 ft/sec 

0.779 

v = -140 7 t e 2lt( - 7í9)p) sin [2 tü(. 2)] = -157.2 ft/sec 
y = Vn 2 tí = V(765.6) 2 - (1572) 7 = 781.6 ft/sec 

From Table A.l, for M« — 0.6269, = 1 -079 

T 0 = 1.079 = 0.079 (519) = 560°R 


ao= JftX, = V04)(m6X56Ó) = 1160 ft/sec 

a 2 = a 0 2 + (V 2 ) = 1.345 x 10 6 - (.2)(781,6) 2 = 1.223 x 10 6 


a= 1106 ft/sec 



a 


781.6 

1106 


0.7067 
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FromTable A.l,forM = 0.6269: ~ = 1.3065, i = 1.079 

P. X 

rn 

For M = 0.7067 = = 1.400, = 1.101 

P T 


p=— -^p«=|—| (1.3065)(1 atra) - (0.933 atm| 
P 0 PU .4. 1 


T T 

T= — — Ttc - i- 

X. X U.101. 


(1.079X519) =|508.6°R 


11.2 The results of Fig. 4.5 are for low-speed, incompressible flow. Henee, from Fig. 4.5, 
at a = 5°, ata = 5°. 


c f =0.75 


c„ = 


0.75 


a/i-mí Vi - (o.6> 2 


10.9381 



11.3 


C P = 


-034 


^mT Vi-(5 8) 3 


—034 

0.8146 


•0.663 


00 C p 


í>o 




+ 




1 + Jl-Ml 


-034 

0.8146 + 

0.3364 " 

(-034) 

1 + 0.8146 

2 


C p = |-0.706 3| 


(c) C p 


p« 




c P - 


-034 


0.8146 + [0.3364(1.067) /1.6292](-034) 
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1-0.7763 


Note the differencesr There is a 17% discrepancy between the three coinpressibility 
corrections. Of the three, experience has shown the Kaiman-Tsien rule to be more accurate. 


11.4 For the pressure coefficient on the airfoil: 



M» 0.3 0.4 0.5 0.6 


0.7 0.8 


C p -0.43 -0.447 


-0.473 -0.513 -0.574 -0.683 
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11.5 When M = M cr , then p at the mínimum pressure poínt is clearly p cr - 


_P_ 

P» 





= (0.528)(1.524)= t0_805l 



Evaluated 
atM -1 


-n 

Evaluated 
at M — 0.8 


11.6 From Appendix A: 


Po 


ForM«> = 0.5, — = 1.186 
P» 


For M = 0.86, ^ = 1.621 


C - P~P« - i*-P- - 2 

q- ^ 


VP„ ) 


= P sllsa. = - 0.7316 

P- P„/P L621 


c ? 


mm 


= (0.7316 — 1)— 1- 1-531 


Check: Using Eq. (11 -58) 


C p = 


y Mi 


^1 + ——-M 2 ^ 
2 


y ^ ^ /fí 


1 + 4- M 

v 2 


T 

r -1 


-1 


(1-4X0# 


( 1 + 02(0# ^ 1 


- 1.53 

U + 02(0.86) 2 J 




It checks! 
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11.7 First, calcúlate C p 0 at point A from the infoimatíon in Figure 11.5(a). The actual 
pressure coefficient is 


rK Ip- ) 


where 


Pa = Pa. Po. 

P- Po P» 

From Appendix A (iuterpolating between entries for more accuracy for this problem). 


ForM m = 0.3: =1.064 

P. 


For M a = 0.435: ^- = 1.139 


Thus, 


Cp y A 


Pa 


f 1.064 


(1.4)(0.3) 2 \ 1.139 
From the Prandtl-Glauert rule. 


-1 =-1.045 


C p , 0 = Cp >A -Jl-K = (—1.045) Vi — (0.3) 2 = - 0.9969 
Forthe case of parí (c) where M» — 0.61, againusing the Prandtl-Glauert rule, 


Opa ~ 


-0.9969 


= - 1.258 


Vi "MÍ Vl“(0-61) 2 

To fínd the local Mach number, Ma, which corresponds to this valué of C P , A , note that 


Cp^A = —11 


Y vp 


or. 
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iX , y - 0-4X0.61)^(-1.258) +1 , 06723 

P. 2 2 

However, 


Pa = Pa Po. where P^ for M» = 0.61 is 1.286 

Paj Po P tet PíO 

Thus, 

j p > /p, = a6723 = 0 j23 

Po Po 1 P- 1286 

Henee, 


^ — 1.912. 
Pa 


From Appendix A, for — — 1.912, Ma — |1 -01 

Pa 


This is cióse enough. Henee, given the numbers in Figure 11.5(a), the numbers in Figure 
11.15(c) are consxstent with the laws of physics. 


11.8 There is a three-dimensional relieving effect for the flow over a sphere. The flow over 
a cylinder is two-dimensional - in order to get out of the way of the cylinder, the flow can 
inove only upwards or downwards. This means it must greatly accelerate to get out.of the 
way of the cylinder. In contrast, the flow over a sphere is three-dimensional — it can move 
not only upward or downward but also sideways. This extra degree of freedom means that 
the flow does not have to speed up so much in flowing over the sphere. Henee, the 
freestream Mach number of the sphere is higher in order to achieve sonic flow on the sphere 
- i.e., the critical Mach number is higher. 
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CHAPTER12 


12.1 Consider a = 5 o = 0.0873 rad. 


4a 


4(0.0873) 


■jMl-1 V(2.6) 2 -l 


0.1455 


From exact theory (Prob. 9.13): c f = 0.148 


0/ 0.148-0.1455 1AA , íno/ 

% error =--x 100 = 1.69% 

0.148 


c d z 


4 a 2 


y/Mi-1 


= c ( g = (0.1455X0-0873) = |0.0127| 


From exact theory (Prob. 9.13): cy = 0.0129 

0.0129-0.0127 , ‘ , - . 

% error =-x 100 = 1.53% 

0.0129 

(b) a = 15° = 0.2618 rad 
4a 


c, = 


= 0.43d 


t/mI-i 

From exact theory (Prob. 9.13): c f = 0.452 

% error = 0452 ~°- 426 x 100 = 3.47% 
0.452 

c d = c e a = (0.436)(0.2618) = 0.114 
From exact theory (Prob. 9.13): Cd = 0.121 


0.121-0.114 1AA 

% error =- 1 -x 100 = 5.7% 

0.121 


(c) a = 30* = 0.5236 rad 
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% 


4a _ 4(0.5236) 

VMÍ-1 V(2-6) 2 -l 


D.873 


From exact theory (Prob. 9.13): c f — 1.19 


% error - 


1.19-0.873 

1.19 


x 100 = 26.7% 


c a = c, oo = (0.873)(0.5236) = |0.457¡ 


From exact theory (Prob. 9.13): c d — 0.687 


% error = 


0.687-0.457 

0.687 


33.5% 


Conclusión: At low a, linear theory is reasonably accurate. However, its accuracy 
deteriorates rapidly at high a. Tbis is no surprise; we do not expect linear theory to hold fór 
large perturbations. It appears that linear theory is reasonable to at least 5 o , and that it is 
acceptable as high as 15°. At 30° it is unacceptable. Keep in mínd that the aboye comments 
pertain to the liít and wave drag coefficíents only. They say nothing about the accuracy of 
the pressure distributions themselves. 


12.2 


2 


(a) C p 


P-P«_ P-P* 


-p M 2 

r oo o 


Y Mi 


í _ ' 




P _ Y C p 


+ 1 


C p = ± 


26 


= + 


26 


: + 


26 


jMl-1 V (2-6) 2 — 1 2.4 


or. 


C p = ± 0.83330 


p __ X M 2 C p + (1.4)(2.6) 2 (0.8333) 




— = ± 3.9430 + 1 


Henee : Examining the physical picture: recalling a = 5 C 


= -3.943 (.0873) + 1 = (0.6558 
P» 


From exacttheory (Prób. 9.13): — = 0.7022 

P. 


0, 0.7022-0.6558 1/vv 

% enror =-x 100 = 6.6% 

0.7022 


+ 3.9430 + 1 = 3.943 (.0873) + 1 = (1.344 
P. 


From exacttheory' (Prob. 9.13): = 1.403 

P« 


1403-1344 

% error - x 100 = 4.2% 

1.403 


= 0.873 rad. 
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(b) For a= 15° = 0-2618 rad: 


Ea = -3.9430 + 1 = -3.943 (.2618) + 1 = -0.0322 (physically impossible) 
P- 


The result from exact theory (Prob. 9.13) is —— 0.315 

P« 


P 3 _ 


3.9430 + 1 = 3.943 (.2618) + 1 = R.032| 


From exact theory (Prob. 9.13): —-— 2.529 

P- 


ft/ 2.529-2.032 1Q70/ 

% error =-—-- x 100 = 19.7% 

2.529 

(c) For a = 30° = 0.5236 rád 


El = -3.9430 + 1 = -3.943 (0.5236) +1 = -1.064 (physically impossible) 
P- 


The result from exact theory (Prob. 9*13) is — ~ 0.072o 

P« 


El. = 3.9430 + 1 = 3.943 (0.5236) + 1 = |3.065l 
P* 


From exact theory (Prob* 9T3): 


■El. 5.687 
P- 


5.687-3.065 

% error = ——y- 

5.687 


x 100 = 46% 


Conclusions: (1) Pressures predicted by linear theory rapidly become inaccurate as a 
mercases. (2) Pressures predicted by linear theory are reasonable only at low valúes of a, 
say belów 5 o . (3) At each valué of a, the % error is much greater for pressure tiran for lift 
and wave drag coefficients. (SeeProb. 12.1). Henee, linear theory works better for c r and 
c a than it does for p. What happens rs that the inaccuracies in p on the top and bottom 
surfaces tend to compénsate, yielding a more accurate aerodynamic forcé coefficient. 


113 



12.3 -E- = — M * C> + 1 wbere C p = ± 2g 

P- 2 -JmITí 


2i9 

Cp = ± =±0.70710 


a /( 3) 2 “1 


P _ ± (1.4)(3) 2 (0.7071)^ [ 1 


— =±4.4550+ 1 



Surface 2: 0 = 5° = 0.08727 rad. 


— = -4.455 (.08727) + 1 = 0.6112 
P» 

Surface 3:0 = 25° = 0.4663 rad 


— = -4.455 (.4363) + 1 = -0.9439 
P- 


Surface 4: 0 = 25° = 0.4363 rad 


— = 4.455 (.4363) + 1 « 2.944 
P- 


Note : Although a negative pressure is 
Dot physically possíble, in order to 
calcúlate the net forcé, we must carry 
ít as süch. 
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Surface 5: 9 = 5° = 0.08727 rad 


-^-=4.455 (.08727)+! =1.3888 
P» 

C[ = —-— — c 0 s25°+í——cos5° (From Prob. 9.14) 

r c vp* p«J vp„ p«J 


c . =---- f(2.944 + 0.9439) eos 25° + (1.3888 - 0.6112) eos 5 o ] 

(1.4)(3) 2 c 

c, = 0.682 —. However 7 — = 0.5077 (From Prob. 9.14) 
c c 

c f =(0-682)(.5077) = 10.346 
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CHAPTER13 




v,= (Mi) = 26.38° 

K_= 9 +v = 20 + 26.38 = 4638° 

Atpoint2: 

a 2 = jftX, = V(¿4)(287)(288) = 340 m/sec 
V 2 - 680 m/sec 
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V 

M?= — 
a ! 


680 

340 


e 2 =o° 

v 2 = 26.38° 
K+^e-v *-26.38° 


At point 3: 

03 = V2 [K_)i + (K + ) 2 ] = V2 (46.38 -26.38) = 10° 
v 3 = V2 [Kj! + (K + ) 2 ] = Vz (46.38 -26.38) = 36.38° 
M 3 = 2.4 

To obtain the other flow variables at point 3, note that: 


^ =7.824 and = 14.62 
Pi P 3 


„ _ P 3 Po, Pe, 

P3-Pl 

Pe, Po, Pl 


14.62. 


(1)(7.824)(1 aten) = jQ.535 aten) 


T T 

—= 1.8 and =2.152 


T 3 = 


t T T 

4-1 X 0 X O 


T 03 T #1 t, 


1 


2.152 


° 5 °' Ti= —— I (1X1.8X288) = |240.9°K 


a 3 = = V(1.4)(287)(240.9) «211.1 m/sec 

V 3 = M 3 a 3 = 2.4 (311.1) = 746.6 m/sec 


u 3 = V 3 eos 0 3 = 746.6 eos 10° = |735.3 m/sec| 


•v 3 = V 3 sin 0 3 = 746.6 sin 10° = [129.6 m/sec 
To lócate point 3: 

Along the C+ Characteristíc: 
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9ave = V 2 (02 + 03 ) = V2 (O + 10) = 5° 

q avc = Vz 0i 2 + JX3) = y 2 (30° + 24.62°) = 27.31° 

— = Tan (0 ave + |i aV e) = Tan (5° + 27.31 °) = 0.6324 
dx 

Thus: 

y = 0.6324 x-0.00765 (1) 

Along the C. characteristic: 

Oave = V2 (0, + 03) = % (20° + 10°) = 15° 

= % (fi, + ni) = % (30 + 24.62) = 27.31° 

— = Tan (9 ave * q ave ) = Tan (15° - 27.31°) = -0.2182 
dx 

y = -0.2182 x + 0.0684 (2) 

Poínt 3 lies atthe intersection of Eqs. (1) and (2) 
y = 0.6324 x-0.00765 
y =-0.2182 x +0.0684 
Solving simultaneously: x = 0.0894 

y = 0.0489 



CHAPTER 14 




Foro. = 30°: 


C p; =2 sin 2 30° = 0.5 


c, = 0.5 eos 30° = p.433 


Cd ~ 0.5 sin 35° = 10.25 


(b) Using modified Newtonian: 


C p — C n sin 2 a 


_ Po Peo P O Pa 




y M 


r p, 

«j v F* co 


(l-4)(2.6 f 


(9.181-1)= 1.729 


For a = 15° 


C =1.729 sin 2 5 o = 0.0131 


c f =0.0131 eos 5° = |0.0131 


c d = 0.0131 sin 5° = [0.0011 41 
For a = 15° 


C 3 =1.729 sin 2 15° = 0.1158 

P.T 


c, =0.1158 eos 15° — k).l 119 


c d = 0.1158 sin 15° = |0.030 
For a = 30° 

C p /= 1.729 sin 2 30° = 0.4323 


c £ =0.4323 eos 30° = ¡0.374) 


c d = 0.4323 sin 30° = p.21S 
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Comparison: 



Exact c f 

Newtonian 


Mod. 

Newtonian 


a 

(Prob. 9.13) 

c, 

% error 

c í 

% error 

5° 

0.148 

0.0151 

90 

0.0131 

91 

15° 

0.452 

0.129 

71 

0.1119 ' 

75.2 

30° 

1.19 

0.433 

63.6 

0.374 

68.6 


Exact c d 

Newtonian 


Mod. 

Newtonian 


a 

(Prob. 9.13) 


% error 

c d 

% error 

5 o 

0.0129 

0.00132 

90 

0.00114 

91 

15° 

0.121 

0.0347 

71 

0.03 

75.2 

30° 

0.687 

0.25 

63.6 

0.216 

68.6 


Conclusión : Newtonian theory gives terrible results for a fíat píate a modérate a at low 
Supersonic Mach numbers. 



From Newtonian theory: 

Cp — 2 sin 2 a = 2 sin 2 20° = 0.234 


121 


c, =0.234 eos a = 0.220 


c ¿ ~ 0.234 sin a = O.0& 


From shock-expansion theory: 

On the ton snrface : V 2 — vj + 9 = 116.2 + 20 = 136.20 

Thís is beyond the máximum expansión angle. Henee, a “void” exists on the top snrface, 
i.e., p 2 = 0. 

On the bottom snrface : From the 9-p-M diagram, 

0 = 24.9° 


M n( = Mi sin p = 20 sin 24.9° = 8.4 


P3 _ 


82.15 


FromProb. 9.13: 


c. 


y M ra Kvt p, 


2 í p 3 P2 


eos a 


and 


sin a 

Cd — o ( 

cosa 


c. =-— (82.15 - 0) eos 20° = 0.2757 

(1.4)(20) 2 

c d = 0.2757 Tan 20° = 0.100 


Forc f : % error ; 


0.2757-0.220 

0.2757 


= 20 % 


For Cd: % error = 


0.100-0.08 

0.10 


= 20 % 


Flote: Newtonian theory works mnch better for blunt bodies, i.e., for large valúes of 9. 
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14.3 


AL -r^O 



(a) Use Eq. (14.7) to estímate the pressure at point A. We ffrst need to obtain 
C pmax , which is a function of p o yp®. From Appendix B for MU = 20, p o yp® = 0.5155 x 
10 S . Henee, 


r 


^1-1 =- 1 — (515.3-1)= 1.837 

V P„ J (l-4)(20) 2 


From Eq. (14.7), at point A on the suiface 


C p = C p , max sin 2 0 = (1.837) sin 2 20° = 0.2149 


Pa 


Since 


Pa 


Y M 2 


/_ \ 


vp« J 


then, 

Pa T MjC t> r t _ (L4X20) 2 (02149) | y ^ 
P. 2 2 


Henee, 


Pa = 61.17 (3.06) = 1187.2 lh/ft 2 

(b) The stagnation temperatüre is found from Eq. (8.40) 


rp -I 

= 1 + tZL Mj = 1 + 0.2 (20) 2 = 81 


Ass umin g an isentropic flow from the stagnation point to point A, 
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or 


Pa / P= 


P 0> 2 vTj 


r -1 


_ Í61.17VT = 
V 515.5/ 


(0.1187) 0-2857 = 0.5439 


T 0 VT» 


T f T 

T a = — I ~ I T« = (0.5439)(81)(500) = j22,028°R| 


(Please note. Relativé to our discussion in Problems 8.17 and 8.18, we know this 
estímate of T A to be too large because we are not taking into account the effect of 
chemically reacting flow.) 

(c) At poínt A. for an isentropic flow, p 0iA = po ,2 


o. A 


y — 1 

1+^-M 

V 2 n J 


r-' = /p, ^ 515.5 

p A /p a 61.17 


= 8.427 


1 + M 2 = (8.427) = (8.427) 0-2857 = 1.8385 


M a 2 = (1.8385-1) -2— = (0.8385)(5) = 4.1925 
Y -1 


M A = 2.05 


(d) a A = “ 704)0716X2X028) = 7275 ft/sec 


V A = a A M a = (7275)(2.05) = [1,49 x 10* ft/sec 


Note : Once again, this estímate of Y A is too high because T A , henee a A , is too high. 

Also note : The purpose of this problem is to illustrate that, from the Newtornan sine- 
squared law for pressure variations, the other-flow field quantities can also be obtained. 
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CHAPTER 15 


15.1 


oO 


¿- / r-¿'_ / _ ¿1 

~ i U-Uf 


oO 


U -o 


“V / / v / / / 




(a) Since the plates are infinite in length, n = u(y) only. Also, v = 0, i.e., the flow is 
in the x-direction only. The goveming equation is Eq. (15.18a), which reduces to the 
following u — u(y), v = 0 and p — const. 


o-A((i—> 

dv dv 


Integrating: 


du 

p. — = const = C} 
dv 


qu « cjy + c 2 
At y = 0, u = 0: c 2 = 0 
Aty = h.u = u e : qu e = c¡h 


h 


Thus: 


/¿u £ 

qu= —- y, or 



The velocity vaxiation is linear between the plates. 


dy h 
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du u. 

t = P — -p — 
dy h 


JlJI] T o + n °_f 320 V 2 288-16 + 110 , 1 Qg/] 
ju a ItJ T+110 V 288.16/ 320 + 110 


p = 1.084 p 0 = 1.084 (1.7894 x 10" 5 ) = 1.94 x 10 


v5 


kg 


m sec 


t = (1.94 x 10 ) 


-S fI _ le o “v .. i 7v3~' 


Vo.oi 


5.82 x 1(T N/m' 


The shear stress is constant. and henee is the same on the top and bottom walls. 


15.2 



u = u(y),v = 0,p = p(x) 


0 = 


¿p + A 

dx dy l 


du^ 
dy J 




du 

dy 



y + ci 


lu,= '(í )^ +c,y+<:2 

f dp^ y 2 c. c, 

u = - _fl + _L y + _L 

Vdx/ 2 /j ¡j ¡i 

At y = 0, u = 0. Thus ci ~ 0 
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At y = h, u - 0. Unís, 




dx7 2 ¡i 


+ A h 

M 



h 

2 


u = - 


dpW 

dxJ 2 fi 






A= ÍA) L + ( A| A 

dy VdxJ ¡i VdxJ 2/i 

On the bottom píate, y — 0: x = p A 

dy 


The velocity profile is parabolic. 


x = r í Al — + f Al A i = — ( dp l 

VdxJ p VdxJ 2/2 ^ 2 VdxJ 


On the top píate, y — h: x — p (——) slnce dy is negative, i.e., the distance away from the top 
píate is in the downward (negative dixection) 


t = p [+ 



h 

A 




^ _ h í dp'j 

T 2 vAJ 

For both the top and bottom walls, 

T-Í(ÍE) 

2 Vdx2 

Shear stress vanes linearly with the magnitude of the pressure gradient. 


Note: Díte to the contení of chapters 16,17, and 18, no hómeworkproblems are 
required. 
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CHAPTER19 


19.1 1 mi/hr = 0.4471 m/sec 


V.» 141 


mi 

h 7 


0.4471 m/sec 


1 mi/hr 


= 63.04 m/sec 


Re c - = (1.23)(63.04X lj) = 6 . 9 . x , 0 4 

M* 1-7894 x 10' 5 

WC-JP-, 1328 ■— 5.04x30^ 

V Re c V6.93 X 10 6 

Noting that drag exists on both the bottom and top surfaces, we have 

■A = 


Dr = 2q TO S C f = 2(V2)(1.23)(63.04) 2 (9.75)(1 -6)(5.04 x 1Q- 4 ) = ]38.4N| 


MO - m . - 0074 ... -3.17X10* 


D,= 


(Cf Xurb 


Re e 1/5 (6.93x10 6 ) 1/5 

3.17 x 10" 3 


(38.4) = 


5.04 x 10" 


(38.4)= [241.5N . 


(Cf-)i am 

Note thatturbnlent skin friction is 6.28 times larger tban the laminar valué. 


192 ( a ) = 3 Q4 x 10^m= |0.304 cm 

V6.93 x 10 6 


= 037x = JU7XL6L , 2 . 54 x 10 - m = PÍ 
' ' Re ' (6.93x10 6 ) 1,5 


era 
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q„ = V 2 (1.23)(63.04) 2 = 2444 N/m 2 


Re c = 5x 10 s 


P*V.(x,-x 0 ) 

P- 


(Xi-X 0 ) = 


5 X Ipy„ (Sx 10 5 )(1.7894 x 10' 5 ) 
P„V„ (1.23X63.04) 


= 0.1154 m 


Laminar drag on (xj - Xo): 


C f = - =L 3 , 28 ■ = 1.878 x 10‘ 3 

X5 x 10 5 


Dr=q» S C f = (2444)(0.1154)(9.75)(1.878 x 10' 3 ) = 5.16 N 
Turbulent drag on (xj - x 0 ): 

0.074 


Cf= 


5 \ 1/5 


D f = 


(5 x 10 5 ) 


f 5.36 x 10 


5.36x10" 3 


-3 A 


17 


1.878 x 10“V 


5.16 = 14.73 N 


■ From Prob. 19.1, the turbulent drag on (X 2 — x 0 ) was 241.5 N. Henee, 


Turbulent drag on (X 2 - xi) = 241.5 — 14.73 = 226.8 N 


Total skin friction drag = [Laminar drag on (xj - x 0 )] + [Turbulent drag on (x 2 - Xi)] 


5.16 + 226.8= 232 N 


19.4 At standard sea level: p«= 0.002377 slug/ft 


T« = 519°R 

a» = J}RT = J(1AX1716)(519) =1117 ft/séc 
V QO = Mo 0 a IO = 4 (1117) = 4468 ft/sec 
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P„V,c _ (.002377X4468X5/12) _ } ] 0 x in7 
U 3-7373 x 10 “ 7 


1.328 1328 

Incompressible Cf= C, - - 7 = = ~¡= ■ ■■ == 

0 7^7 VI-18 x 10 7 


= 3.866 x 1o -4 


FromFig. 18.8: 

C/C f * 0.85; Cf= 3.286 x 10' 4 

T o 

D f = q„ S Cf = (y2)(.002377)(4468) 2 (5/12)(3.286 x 10* 4 ) 


Df = 3.248 Ib| on one si de of the píate. 


19.5 For incompressible flow: 

= 0-074 ^-0 : 074-=2.85 x 10 " 3 

” Re c 1/5 (1.18 x 10 7 ) l/5 

FromFig. 19.1: Cf =» 1 . 6 x 10 * 3 

(Tire effect of Mach number is to reduce Cf by about 44% in tbis case.) 

FromProb. 19.4, the laminar valué of D f is 3.248 for a valué of C r = 3.286 x 10 -4 . 
Henee, the turbulent valué is 


D f = 


i" 2.85 x 10~ 3 
V3286 x 10" 


(3.248) = [28.2 lb| 


19.6 FromEq. (18.32): 


¿h ¿h d ( ai 

pu ■— + pv — * — u — 

¿k cfy dj V ¿y7 


From Eq. (18.41) vdth Pr = 1: 


da di 0 d 

pll __¡L +pv — = — 

oí ¿y ¿y 


di 


Eqs. (1) and (2) are identical. Henee 


( 1 ) 


( 2 ) 
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ho = c 5 + c 2 u, where cj and c 2 are constante. 
Atthe wall, u = 0 and h 0 = h = h w . Henee, 

hw - ci + 0 , or ci = hw 
At the boundary layer edge: 

h 0t = C! + c 2 u e = h w + c 2 u,. 



Thus: 


b<, = C] + c 2 u = h w + 



u 


Sínce 


h = c p T, then 



19.7 FromEq. (18.70), 


0.763 Pr 0 65 (p e p c ) 



(baw hw) 


( 1 ) 


where, from Eq. (18.82), the velocíty gradient is giyen by 


_ 1 2(p e -p,) 

dx R Al p■ 


( 2 ) 


The subscript e denotes properties at the outer edge of the stagnation point boundary layer, 
i.e_, p c and p e are the inviscid stagnation point valúes of pressure and density. The speed of 
sound iñ the ambient atmosphere is 


a*- = V(1-4)(287)(246.1) = 314.5 m/sec 

(a) For Veo — 1500 m/sec, we have 
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1500 

3145 


= 4.77 


From Appendix B (nearest entry). 


=29.52 

P« 


and from Appendix A (nearest entry). 


T 

- £ - = 5.512 
T 

rrt 


Henee, 


P 0 .2 = p e = (29.52)(583.59) = 1.723 x 10 4 N/m 2 
T 0 = T e = (5.512X246.1) = 1357 K 


p B = 


Pe 

RT C 


1.723 x 10 4 
(287X1357) 


= 0.044 kg/m 3 


From Southerland’s law, Eq. (15.3), using the standard sea level valué of p 0 
kg/(m)(sec) at T 0 = 288K, we have 


Ae 


T.) M T,+11Q 

tJ T e +110 


(1357' 


V 288 


288+110) 
1357 + 110/ 


= 2.77 


Pe = (2.77)0 -789 x 10’ 5 ) = 4.957 x 10 ' 5 kg/(m)(sec) 

From Eq. (2) above 

dn^ = ! | 2(p e -pJ 1 1 2(1.723 x 10 4 -583.59) _ g ^ 

dx R)( p e (0.0254) V 0.044 


Assuming a recovery factor r = 1, then h a w — h 0 . 

V 2 v 2 C1500) 2 

baw = ho = h<o + —r— = c p Ta, + —f— = (1008)(246.1) + ^ 

2 2 2 

= 2.48 x 10 5 +11.25 x 10 5 = 13.73 x 10 5 joule/kg 


haw = c p Tw = (1008)(400) = 4.032 x 10 5 joule/kg 


1.7894 x 10 ‘ 5 


c 10 4 /sec 
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The ‘ rho-mu” product is 


Pe^ = (1,044)(4.957 x J O' 5 ) = 2.18 x 10 


.-6 


(kg) 2 


m sec 


From Eq. ( 1 ) above 


q w =0,763 Pr " 0,65 (p^) (h aw ^h w ) 

ax 


0.763 (0.72)‘° 65 (2.18 x 10" 6 )(3.42 x 10 4 ) % (i: 
joules 


= 369.3 


sec(m 2 ) 


569.3 


watt 


m 


(b) For Y*, = 4500 m/sec, we have 


^=— = ■^• = 14 . 3 ! 
a„ 314.5 


From Appendíx B (interpolated) 


— =264.0 
P» 

From Appendíx A (interpolated) 
T 

— = 41.94 
T 

eo 


Thus: 


p c = (264)(583.59) = 1.54 x 10 5 N/m 2 
T e = 41.94 (246.1) = 10,321 K 
1.54 x 10 J 


Pe = 


RT c (187)(10,321) 


- = 0.052 kg/m" 


_í TV' 2 T 0 + 11Q _(' 1032lV ’ /2 r 288 + 110 ) 
Mo Vlj T c +110 V 288 J U0321 + H07 


.73 - 4.032) x 10 s 


= 8.186 
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p e = (8.186X1.7894 x 10' 5 ) = 1.465 x 10^ kg/(m)(sec) 

H; 

From Eq. ( 2 ) 


ÍL-i Í 2 (P-"P-) —L- jjflÜ xjo»-5Bao , 956xloVc 

dx R| (0.0254) V 0.052 

h aw = !»«,+ -^—=2.48 x 10 5 + = 1.037x 10 7 joules/kg 

PeUe - (0.052)(1 -465 x 10- 4 ) = 7.62 x 10' 6 

m sec 


q w =0.763 Pr ^ 5 (p e p e ) , —2- (h aw -h w ) 

dx 


du 


= 0.763 (0.72y 0 65 (7.62 x 10‘ 6 )(9.56 x 10 4 ) % (1.037 x 10 7 - 4.032 x 10 5 ) 


2.218x10 


4 watts 


m 


Comparing the resnlts from parts (a) and (b), we note 



2218 x 10 4 
369.3 


= 60 


When the velocity increased from 1500 m/sec to 4500 m/sec, a factor of 3, the heat transfer 
increased by a factor of 60. This illustrates the rapid growth of the importance of 
aerodynamic heating as vehicles fly faster. well into the hypersonic flight regime. A simple, 
approximate analysis for aerodynamic heating which assumes very high Mach numbers (so 
that h aw » h w ) indícates that aerodynamic heating is proportional to V, 3 . (See for example, 
Anderson. Introdnction of Flight, 4* ed., McGraw-Hill, 2000, page 570.) For the present 
example, in going from a relatively low, not quite hypersonic condition (M„ = 4.77) to a 
relatively high Mach number of M® = 14.31, the increase was eveh faster. 
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